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Abstract 

We consider weak solutions to a two-dimensional simplified Ericksen-Leslie system of compressible 
flow of nematic liquid crystals. An initial-boundary value problem is first studied in a bounded 
domain. By developing new techniques and estimates to overcome the difficulties induced by 
the supercritical nonlinearity |Vdpd in the equations of angular momentum on the direction 
field, and adapting the standard three-level approximation scheme and the weak convergence 
arguments for the compressible Navier-Stokes equations, we establish the global existence of 
weak solutions under a restriction imposed on the initial energy including the case of small 
initial energy. Then the Cauchy problem with large initial data is investigated, and we prove the 
global existence of large weak solutions by using the domain expansion technique and the rigidity 
theorem, provided that the second component of initial data of the direction field satisfies some 
geometric angle condition. 

Keywords: Liquid crystals, compressible flows, weak solutions, Galerkin method, weak 
convergence arguments. 
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1. Introduction 

We consider the existence of weak solutions to an initial-boundary value problem for the 
following two-dimensional simplified version of the Ericksen-Leslie model in a bounded domain 
n C which describes the motion of a compressible flow of nematic liquid crystals: 

9ip + div(pv) = 0, (LI) 
dt{pv) + div(pv ® v) + VP(p) = pAv + (p + A)Vdivv 

-z/div(^Vd0 Vd- i|Vd|2lj, (L2) 

(9jd + V- Vd = e(Ad+ iVdl^d), (L3) 

where p is the density of the nematic liquid crystals, v the velocity and -P(p) the pressure, 
d G := {d G I |d| = 1} represents the macroscopic average of the nematic liquid crystal 
orientation field. The constants p. A, u, and 9 denote the shear viscosity, the bulk viscosity, the 
competition between kinetic energy and potential energy, and the microscopic elastic relation 
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time for the molecular orientation field, respectively. I denotes the 2x2 identical matrix. The 
term Vd© Vd denotes the 2x2 matrix whose {i, j)-th entry is given by dx^d-d^jd, for 1 < i, j < 2, 
i.e., Vd © Vd = (Vd)'''Vd, where (Vd)''' denotes the transpose of the 2x2 matrix Vd. 

Nematic liquid crystals are aggregates of molecules which possess the same orientational 
order and are made of elongated, rod-like molecules. The continuum theory of liquid crystals 



was developed by Ericksen [5[ and Leslie [21| during the period of 1958 through 1968, see also 
the book by de Gennes [9|. Since then there have been remarkable research developments in 
liquid crystals from both theoretical and applied aspects. When a fluid containing nematic liquid 
crystal materials is at rest, we have the well-known Ossen- Frank theory for static nematic liquid 



crystals. The reader can refer to the work by Hardt-Kinderlehrer-Lin [llj on the analysis of 
energy minimal configurations of nematic liquid crystals. In general, fluid motion always takes 
place. The so-called Ericksen-Leslie system is a macroscopic continuum description of the time 
evolution of materials under the influence of both flow velocity field v and macroscopic description 
of the microscopic orientation configurations d of rod-like liquid crystals. 

In 1989, Lin 25| first derived a simplified Ericksen-Leslie system modeling liquid crystal flows 
when the fluid is incompressible and viscous. Subsequently, Lin and Liu |27|, l28| established some 
analysis results on the simplified Ericksen-Leslie system, such as the existence of weak and strong 
solutions and the partial regularity of suitable solutions, under the assumption that the liquid 
crystal director field is of varying length by Leslie's terminology, or variable degree of orientation 
by Ericksen's terminology. 

When the fluid is allowed to be compressible, the Ericksen-Leslie system becomes more com- 
plicated. From the viewpoint of partial differential equations, the system fll.ip - fll.3p is a highly 
nonlinear hyperbolic-parabolic coupled system, and it is challenging to analyze such a system, 
in particular, as the density function p may vanish. It should be noted that the system (11.11) - 
(11.31) includes two important subsystems. When p is constant and v = 0, the system (ll.ip - (ll.3l) 
reduces to the equations for heat flow of harmonic maps into on which there have been exten- 



sive studies in the past few decades (see, for example, the monograph by Lin and Wang [29j and 
the references therein). When d is constant, the system (11.11) . (II. 2p becomes the compressible 
Naiver-Stokes equations, which have attracted great interest in the past decades and there have 
been many important developments (see, e.g., for a survey of recent developments). 

Since the supercritical nonlinearity |Vdpd causes significant mathematical difficulties, Lin in 



25l | introduced a Ginzburg-Landau approximation of the simplified Ericksen-Leslie system, i.e.. 



|Vdpd in (II. 3p is replaced by the Ginzburg-Landau penalty function (1 — |dp)/e or by a more 
general penalty function. Consequently, by establishing some estimates to deal with the direction 
field and its coupling/interaction with the fiuid variables, a number of results on the Navier-Stokes 
equations can be successfully generalized to such Ginzburg-Landau approximation model. For 



examples, when p is a constant, i.e., the homogeneous incompressible case, Lin and Liu [27 
proved the global existence of weak solutions in two and three dimensions. In particular, they 
also obtained the existence and uniqueness of global classical solutions either in two dimensions 
or in three dimensions for large fiuid viscosity p. In addition, the existence of weak solutions 
to the density-dependent incompressible fiow of liquid crystals was proved in [l3, [s^. Recently, 
Wang and Yu [38j , and Liu and Qin 33| independently established the global existence of weak 
solutions to the three-dimensional compressible fiow of liquid crystals with the Ginzburg-Landau 
penalty function. 

In the past a few years, progress has also been made on the analysis of the model (11.10 - 
(11.30 by overcoming the difficulty induced by the supercritical nonlinearity |Vdpd. For the 



incompressible case, the existence of weak solutions in two dimensions was established in [26 
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and the local and global existence of small strong solutions in three dimensions was proved in 



12|, |23|, l30|, |37| . For the compressible case, the existence of strong solutions have been investigated 



extensively. For examples, the local existence of strong solutions and a blow-up criterion were 



obtained in [15|, |16|, while the existence and uniqueness of global strong solutions to the Cauchy 



problem in critical Besov spaces were proved in [IJ] provided that the initial data are close to an 



equilibrium state, and the global existence of classical solutions to the Cauchy problem was shown 
in 



22l | with smooth initial data that has small energy but possibly large oscillations with possible 
vacuum and constant state as far-field condition. To our best knowledge, however, there are no 
results available on weak solutions of the multi-dimensional problem f ll.ip -( lL3|) with large initial 
data due to the difficulties induced by the compressibility and the supercritical nonlinearity. It 
seems that the only global existence of large weak solutions to fll.ll) - fll.3p was shown in the 
one-dimensional case in j^. 

In this paper, we will establish the global existence of weak solutions to the two-dimensional 
problem fll.ip - fll.Sp in a domain C with initial conditions: 

p(x,0) =po(x), d(x,0) = do(x), (pv)(x,0) = mo(x) in Q, (1.4) 

and boundary conditions: 



d(x, t) = do(x), v(x, t) = 0, x G dQ, t > 0, if is a bounded domain; 
(p, V, d)(x, t) — 7- (poo, 0, 62) as |x| — 7- 00, if = M^, 



;i.5) 



where Poo > is a constant and 62 = (0, 1) G is the unit vector. For simplicity, we call fll.ip - 
fll.51) the initial-boundary value problem when i7 G is a bounded domain, and the Cauchy 
problem in the case of f2 = M^. 

Before stating our main result, we remark that the constants p. A, u, and 6 satisfy the physical 
conditions: 

p>0, A + p>0, > 0, 9>0. 
And the pressure -P(p) is usually determined through the equations of state, here we focus our 



study on the case of isentropic flows as in [38| and assume that 

P(p) = Ap^, with A>0, 7 > 1. 
In addition, for the sake of simplicity, we define 

/:=/t:=(0,T), QT = nxI, 

jr(i) := jr(p^v,d) := / Vv^ + (A + p) |divv|2 + ^(| Ad + | Vd^dp)) dx, (1.6) 

Jn 



and 



f /l Imp z/^lVdPA 
g(t):=g(p,m,d) = y^(-^l|^>o} + g(p)+ ' Jdx with m = pv, (1.7) 

where l{p>o} denotes the character function, and 

Q{s) = (s^ — 'ysp]^^ + (7 — 1)pL) some constant poo > 0. (1.8) 

7 — 1 
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Obviously, Q{s) > for any s > 0, and £{t) > 0. Throughout this paper, we use the bold fonts 
to denote the product spaces, for examples, 

LP{n) := H^(fi) := {H^{n))^ = (Vro''(^))', H'^(^) ■= (W^'^n))^; 

and the Sobolev space with weak topology is defined as 

C'il K^^^m := |f : / ^ ^ f • gdx G C{I) for any g G L^(f])| . 

Our first result on the existence of weak solutions in a bounded domain Q reads as follows. 

Theorem 1.1. Let the constant poo > 0, f2 6e a hounded domain of class C"^'" with a G (0, 1), 
and the initial data Pq, iuq, do satisfy the following conditions: 

Q{po)eL\Q), po > a.e. inQ, (1.9) 

I 1 2 

27 nig 1 

mo G LT+l(^]), mol{p„=o} = a.e.inQ, l{po>o} G L (fi), (1.10) 

Po 

|do| = 1 inQ, do(x) G H2(^]). (1.11) 

Then, there exists a constant Cq, such that if £q := £{Q) satisfies 

£o < Cov, (1.12) 

the initial-boundary value problem M.1\) - [T7B^) has a global weak solution {p,v, d) on I = It for 
any given T > 0, with the following properties: 

(1) Regularity: 

0<p a.e.mQr, p e C\l Ll^,^{n) n C\l L^n)) n ^-^^ {Qt) , (1.13) 
V G L2(/,HJ((])), pv G L-(/,L^(n)) nC°(J,L2:,(n)), (1.14) 

|d| = 1 a.e. zn Qt, d e {1 ,11^ {n)) r} C\I {n)) , dtd e L^I ,V {p)) , (1.15) 
d\dn = doloQ G C^'^{dri) in the sense of trace for a.e. t E I, (1-16) 

where p E [1,7), and rj G (0,7— 1) . 

(2) Equations U-l\) and U.^) hold in {V'{Qt)Y, and equation holds a.e. in Qt. 

(3) Equation / li.ij) is satisfied in the sense of renormalized solutions, that is, p,v satisfy 

dth{p) + div[6(p)v] + [ph\p) - h{p)] divv = m V\n^ x J), 
provided fp,vj zs prolonged to be zero on \ Vt, for any h satisfying 

b G C°[0, 00) n C\0, 00), \b'{s)\ < cs-^«, s G (0, 1], Ao < 1, (1.17) 
and the growth conditions at infinity: 

21 — 1 

\b'{s)\ < cs^\ s > 1, where oO, 0<1 + Ai< —- — . (1.18) 
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(4) The following finite and bounded energy inequalities hold: 



^^ + T{t)<0 inV'il), (1.19) 



£{t) + f J^(s)ds < 8q for a.e. te I. 
Jo 



(1.20) 

Remark 1.1. It should be noted that for the incompressible case the constant Cnin a small 
energy condition similar to (I1.12p depends on the domain fi, since the authors in used the 
following interpolation inequality: Given a bounded domain i7 C M^, there exists an interpolating 
constant ci{Q) depending on Q, such that 



||Vd||^4(^) < c^im\^<i\\lin) + II V^dH^.^^^H Vd| 



Inn)) (1-21) 



for any Vd G H^(r2) (see |J, Lemma 2.4]). In this paper, we use another version of interpolation 
inequality (see fl2.12p ). where the interpolating constant takes the value of 2. Thus our constant 
Co in (11.121) is independent of fl. In particular, from our computations (see Proposition 12.11) we 
can take Cq = 1/4096 in (I1.12p . Of course, this value is not optimal. 

Remark 1.2. In the above theorem, we assume that dn G H^(f2) for simplicity. If this condition 
is replaced by do G H^(f2) and do G C^''^(9f2) as in 26|, then the above theorem still holds. In 



addition, the proof of Theorem 11.11 remains basically unchanged if the motion of the fluid is 
driven by a small bounded external force, i.e., when (11.20 contains an additional term pf(x, t) 
with f a small, bounded and measurable function. 

We now describe the main idea of the proof of Theorem II. 1[ For the Ginzburg-Landau ap- 
proximation model to fll.ip - (ll.3l) . based on some new estimates to deal with the direction field 



and its coupling/interaction with the fluid variables, Wang and Yu in [38| adopted a classical 
three-level approximation scheme which consists of the Faedo-Galerkin approximation, artificial 
viscosity, artificial pressure, and the celebrated weak continuity of the effective viscous flux to 
overcome the difficulty of possible large oscillations of the density, and established the existence 
of weak solutions. These techniques were developed by Lions, Feireisl, et al for the compress- 
ible Navier-Stokes equations in [8|, [l9|, l3l| , and we refer to the monograph 3^ for more details. 
Compared with the Ginzburg-Landau approximation model in 38(], however, the system (11.10 - 
(11.30 is much more difficult to deal with due to the supercritical nonlinearity |Vdpd in (II. 3p . 
Consequently, one can not deduce sufficiently strong estimates on d from the basic energy in- 



equality (ll.20p in the general case, such as V^d G L'^{I, L^(r2)) as in [38|; and on the other hand, 
one can not establish the global existence of large solutions to the equation (II. 3p for any given 
V G C"^(J, Cq{D.)) as in the Ginzburg-Landau approximation model. Recently, Ding and Wen 
obtained the global existence and uniqueness of strong solutions to the two-dimensional density- 
dependent incompressible model with small initial energy and positive initial density away from 
zero. In [3], they can deduce V^d G L^{I, L^(f2)) from the basic energy inequality in two dimen- 
sions under the condition of small initial energy. In fact, they first got z/|| Vd||^2(Q)/2 < Sq and 
1/6*11 Ad 11^2 (Q < ^^0 + '^^ll I Vd| |^4^Q^^ from the basic energy inequality. Then, by interpolation in- 
equality (11.210 and elliptic estimates, they immediately obtain || V^d||L2(Qj,) < C{So, || V^do||L2(n)) 
for some constant C{£o, || V^do||L2(n)) depending on the arguments, if initial energy is sufficiently 
small. Motivated by this study, for our compressible case we impose the restriction (I1.12p on the 
initial energy, and consequently deduce the desired energy estimates on d from the basic energy 
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inequality as in [38|. With these estimates in hand, we also find that the three- level approxima- 
tion scheme can be applied to the initial-boundary value problem fll.ip - fll.Sp . if we can construct 
solutions to the following third approximate problem: 

pt + div(pv) = eAp, (1.22) 

/ (pv)(t) • *dx - / mo ■ *dx 

JQ. Jn 

pAv + (/i + A) Vdivv - AVp^ - SVp^ - e(Vp ■ Vv 

IVdpiM 

*dxds 



Jn 



1.23) 



- div(pv (g) v) - z/div |^Vd Vd - 
for all t G / and any functions ^ G X„, 



(9td + v-Vd = 0(Ad+|Vd|M), (1.24) 

where the n-dimensional Euclidean space X„ will be introduced in Section [31 and e,6, (3 > are 
constants. Fortunately, we can show the local existence of strong solutions d G C°(/, H^(f2)) to 
the equation (11.231) with smooth initial data for any given v G {v G C°(/, H^(i7) nL°^(i7)) | dt^v G 
L^(J, H^(n))}, with the aid of the well-known Ehrling-Nirenberg-Gagliardo interpolation inequal- 
ity (see Lemma [3.11) . With the result of local existence in hand, adapting some techniques in 
[34 . 38| . and making use of equivalence of norms in X„ and the operator forms of (v, S^v) (see 



(14.91) and (14.111) for the operator forms deduced from approximate momentum equations (I1.23P ). 
and establishing some additional auxiliary estimates on dt^v, we can apply a fixed point theorem 
in some bounded subspace of the Banach space {v G C([0,T*],X„) | ^t^r G L^((0, T*), H^(i7))} 
for some T* G (0, T] to obtain the local existence of the third approximate problem. Noting that 
the regularity of dtv depends on dtp, this is why we need the properties on dtp in Proportion 13.11 
on the solvability of (ll.22p . Then, we can repeat the fixed point argument to extend the local 
solution to the whole time interval / by establishing the uniform-in-time energy estimates for the 
third approximate problem. 

Next, we state our second result on the existence of global weak solutions to the Cauchy 
problem in with large initial data: 

Theorem 1.2. Assume that, for some positive constants poo and (1^2! (Po? mo) satisfies U.9\) . 
U.10\) with in place of Q, and 

Idol = 1 a.e. m R^ do(x) - eg G H1(M2), ^02 > ^02. (1-25) 

where do2 denotes the second component of do. Then, the Cauchy problem U.l\ )- in~S\) has a 
global bounded energy weak solution (p, v, d) for any given T > with the following properties: 

(1) Regularity: 

0<pa.e. mR^x/, Q(p) G L~(/, L^(R2)), 

p e c°(/, Ll^^n')) n c\i, L^m n L7+''(R2 X /), 

V G L\I, {D^'^iR'))') n L\I, LLm, pv G C'{I, l£,(1^')), 

p|vpGL-(/,Li(M2))nL2(J,LL(K')), 

dtd G Lt(/,LL(M2)), |d| = 1 a.e. m x /, 

d G c°(/, Hi„(R2)), Vd G l2(/, u^r") n l^(r2 X J), ^ ■ ^ 

where p G [1,7), G (0, 7 — 1), g > 1, and Q' C R^ is an arbitrary bounded subdomain. 
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(2) Equations U.l]) . holds in V'(R'^ x I), equation U.3\) holds a.e. m x J. Moreover, 
equation U.l\) is satisfied in the sense of renormalized solutions. 

(3) The solution satisfies the energy equality U.20\) with in place ofQ. 

Remark 1.3. The notation D^'^(M^) in the above theorem denotes the homogeneous Sobolev 
space on M^, i.e., D^''^(U.'^) = C(f (M^)" "^^("^^^ where the symbol "overhne with norm" means 



completion with respect to that norm. Please refer to [3J, Section 1.3.6] for more basic conclusions 
concerning the homogeneous Sobolev spaces. 



Remark 1.4. In Theorem II. 2[ the choice of e2 is for convenience only. In general, one can 
choose any reference vector e G S-*- and require that the image of d is contained in a hemisphere 
around e. Let Poo = and suppose that in addition to the assumptions above we have also 
Po G L^(]R^). Then the result above remains valid, moreover there holds J^2 pdx = J^2 Podx in 
/. In addition, in view of the proof of Theorem 11.21 and the following interpolation inequality: 
there exists a constant c such that 



i^Ba) ^ 4v\\Hi(^Bn)Mh{Ba) any v G H\Br) with R>1, 

the above theorem still holds if we impose a small initial energy condition similar to f ll.l2p to 
replace the geometric angle condition > ^02 f !1.25p . Here 5^ C M? is the ball of radius R. 



Remark 1.5. Recently, Lin-Lin- Wang [26| established the existence of large weak solutions for 
the corresponding incompressible case in a two-dimensional bounded domain. However, it is not 
clear whether the method of localized small energy can be applied to the compressible model, 
due to lack of regularity of ||v||j^oo(-/ L2(f^)). Hence it is still an interesting open problem whether 
a general existence result of large weak solutions holds in bounded domains for the compressible 
case. In particular, we can not generalize the above result of the Cauchy problem to the bounded 
domain case, since it is not clear whether the rigidity theorem (see Proposition 16. 2p holds in the 
bounded domain case. 

Remark 1.6. When we completed this paper, we were informed that Wu and Tan [39*] just 
finished the proof of the existence of global weak solutions to fll.ip - fll.Sp in by using Suen 
and Hoff's method [SSj, if the initial energy is sufficiently small, the coefficients p and A satisfy 
< A+p < ^i|^p, and the initial data (vq, vq) satisfies ||vo||Lp(Ri') + ||do||Lp{R3) <oowithp>6. 
However, it is not clear whether their proof can be applied to (ll.ip - fll.3p in a three-dimensional 
bounded domain. 



Here we also describe the main idea of the proof of Theorem 11.21 Recently Lei, Li and 
Zhang [20] proved the rigidity theorem in (see Proposition 16. 2p and obtained the estimate 
on II V^d||i2(/ L2(R2)) from the energy inequality, if the second component of initial data of d 
satisfies some geometric angle condition. Motivated by this study, we impose the geometric 
angle condition dQ2 > ^02 fll.25p in place of (11.121) . We first establish the local solvability 
of the Cauchy problem (ll.24p on d (see Proposition 16.11) . which can be shown by following 
the proof of the local solvability of d for the problem (ll.24p defined in a bounded domain and 
using domain expansion technique. Then, using the rigidity theorem, elliptic estimates and 
interpolation inequalities in R^, we can follow the proof of global existence of solutions to the 
third approximate problem fll.22p - fl.24p defined in a bounded domain to establish the global 
existence of solutions to the approximate problem (ll.22p . (11.230 defined in a ball Bji and (11.240 
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defined in (see Proposition 16 .41) . Finally, we adapt the proof in SJ, Section 7.11] on the Navier- 
Stokes equations to prove Theorem II. 2[ by using Proposition 16.41 with the mass and momentum 
equations defined in the bounded invading domain Bji and letting i? — >■ oo. 

The rest of paper is organized as follows. In Section [21 we deduce the basic energy equalities 
from fll.ip - fll.3l) and derive more energy estimates on d under the assumption (11.121) . In Section 
[31 we list some preliminary results on solvability of two sub-systems (11.221) and (11.231) in the 
third approximate problem, while in Sections [H we establish the unique solvability of the third 
approximate problem. In Sections [5], exploiting the standard three-level approximation scheme, 
we complete the proof of Theorem II. 1[ Finally we provide the proof of Theorem 11.21 in Section 

m 

2. Energy estimates 

This section is devoted to formally deriving the energy estimates from (ll.ll) - (ll.3p in a bounded 
domain O, which will play a crucial role in the proof of existence. Some of our results have been 
established in 0, [i^] for the incompressible case. 

2.1. Energy equalities 

We consider a classical solution (p, v, d) of the problem (ll.ll) - (ll.3p with initial and boundary 
conditions (ll.4p and (11.51) . First we verify that 



|d| = 1 in Qt, if Idol = 1 in (2.1) 
Multiplying the d-system (11.30 by d, we obtain 

-dt\d\^ + -V • V|dp = 9{Ad ■ d + iVdHdH. 



Since 

it follows that 



2 

A|d|2 = 2|Vdp + 2Ad-d, 



a^ddl" - 1) - eA{\d\' - 1) + V ■ V(|d|' - 1) - 20|Vd|'(|d|' - 1) = 0. (2.2) 

Multiplying (12. 2p by (|dp — 1) and then integrating over fl, we use the boundary conditions to 
get 

- [ (|d|2-l)2dx< [ (4e|Vd|2 + divv)(|d|2-l)2dx 
dt Jn Jn 



<||4^|Vd|^ + |divv||Uoo(f,) / (|d|^-l)^dx. 

Jn 

We assume that (v, d) satisfies the following regularity: 

114611 Vdp + |divv|||ii(j^ioc(Q)) < oo, 

then, using Gronwall's inequality, we immediately verify (12. ip . We shall see that all the couples 
(v„, d^,) in the third approximate solutions constructed in Section [H satisfy the regularity above. 
Similarly, we can verify 

^2 > c?02 if do2 > for some given constant rfo2! (2-3) 
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< 



iVdr + -divv 



L°°(n) 



< 



6'|VdP + -divv 



L°°(n) 



1 \ 

+ 2^02^11 1 Vdllli^(f^) I ||u;~||i2(f^) + -\\Vuj- 



which yields 



— Il^^ Whin) < I 2 



0|Vdp + -divv 



+ ^02^ll|Vd"i2 



L°°(n) 



Hence, by Gronwall's inequahty, one obtains 



(2.4) 



which will play a crucial role in the existence of large weak solutions to the Cauchy problem. 
Here we have denoted the second component of d and do by ^2 and do2, respectively. 
Next, we give the deduction of (12.31) for the reader's convenience. Let 

u = d2 — dQ2, co~ = minju;, 0}, 

we can deduce from (11.31) that 

dtu - eAu = 6*1 Vdp(a; + ^02) - v ■ Vw. 
Multiplying (12. 4p by and integrating over Q, we get 

1 Cj? II — 1|2 

2dt lli2(iR2) + 6^11 11^2(5^) 

= / [6'|Vdp(a; + ^02) - V ■ Vwjw^dx 
Jn 

'^"IIl2{q) + ^02^11 llLoo(f^) / iVw^lcj^dx 

Jn 



ltd 



L2(^Q) < \\UJ (0)||^2(!^) 



2 Jo(2riVd|2 + idivv|| ^^^,„,+dhm^d\\\U^a))^s 



0, 



which implies (12. 3p . 

Now we derive the energy equality. With the help of (12. ip . we can deduce the basic energy 
equality. To this end, we multiply equation (II. 2p by v and integrate over fl to deduce that 



(2.5) 



/ (Vd)' Ad- vdx, 
'n 



where we have used integration by parts, the boundary condition of v, the mass equation (II. ip 
and the equahty 

div (vd Vd - ^iVdpI^ = (Vd)^Ad := {didj)2x2Ad. 
On the other hand, multiplying (II. 3p by —(Ad + | Vdpd) and integrating over Q, we obtain 



- / ^td -Addx- / (v • Vd) • Addx = / |Ad+ |Vd|M|Mx, 
'n Jn Jn 



(2.6) 
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where we have used the fact that |d| = 1 to get 

{dtd + V ■ Vd) ■ I Vd|M = ^(1 Vdp9t|dp + \Vd\\ ■ V\d\^) = 0. 
Using the boundary condition of d, and integrating by parts, we deduce from fl2.6p that 



1 d 
2dt 



Vd|Mx + / |Ad+ |Vd|M|Mx= / (v • Vd) • Addx. (2.7) 
Q Jq Jn 



Consequently, fl2.5p + f l2T|) x;/ gives the energy equahty in the differential form: 

±g{t)+J^{t)=Q, (2.8) 

where the energy £{t) and dissipation term J^(t) are defined by (11.71) and (11.61) . Finally, integrat- 
ing (12. 8 p with respect to time over (0,t), we get the energy equality 

S{t)+ [ J^{s)ds = So, tel, (2.9) 
Jo 

where Sq = S{0) denotes the initial energy. Due to the weak lower semicontinuity of norms on 
the right-hand side of (12. 9p . for weak solutions one expects an inequality rather than an equality. 

2.2. More a priori estimates on d under the condition U.12\} 

In this subsection, we will deduce more estimates of d under the assumption (I1.12p , including 
the case of small initial energy. For simplicity in the deduction, we always use the positive 
constant Ci{fl) to denote various constants depending on Q. 

Multiplying (II. 3p by Ad, integrating the resulting equation over Q, and using integration by 
parts, we find that 

1 d 
2dt 



Vd|Mx + ^ / |Ad|Mx = ^ / |Vd|Mx+ / (v ■ Vd) • Addx. (2.10) 

Q Ju Jci Jci 



Thus, fl23|l + fl230P x;y gives 

-£{t)+ [ [/i|Vv|2+ (A + /i)|divvp + z/^|Adp] dx = z/^ [ |Vd|Mx. (2.11) 
^ Jq Jn 



d_ 

Tt 



Recalling VL C M^, exploiting 36|, Lemma 3.3], i.e., 

Wvfuin) < 2||Vi;||^2(o)||f ||i2(n) for any v G HI{VL), (2.12) 

we can infer that 

ll^jC?illL4(Q) < {\\didQj\\Li(Q) + 2^/%diV{dj - dQj)\\^^2(Q^^\\di{dj - c?oi)||^i^(j^)^ 

< 8 (^||<9irfoj|lL4(j^) + 2\\diV{dj - dQj)\\ljif^^^\\di{dj - (ioi)|lL2(n)^ 

< 32(^ci((])||aicioil|ii(n)(||ai(ioj||^i(f^) + \\didj\\l2^^)) 
+ \\^didj\\l2(^)\\di{dj - c?oi)||i2(f^)) for 1 < 2, j < 2. 
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(2.13) 



Consequently, making use of fl2.9p . fl2.13p . Holder's inequality, and the elliptic estimate 

||V2rf,||L2(n) < 2\\Ad,\\mn) + Q\\V^dih2(^n) (2.14) 
deduced from [lo!, Collary 9.10], we see that the last term in (12. lip can be estimated as follows: 

p9 [ |Vd|Mx<4z/^ V [ l^irfjfdx 



+ \\'^dd,\\Un)\md,-doMHn)) 



(2.15) 



(2.16) 



<128ue (ci(^)i|5.^o,il^H^^)(ll^^^o,lll,i(n) + l|5.rf.lli^(n)) 

1<«J<2 

+ Q\\^'doj\\l^a)\\di{dj - doj)\\l.^n) + 2\\Ad,\\l,^^^\\di{d, - doj)\\h^n) 
<128u9ic^in) + 12)||do||^.(^) (l|do||^2(f,) + ^) + 1024^^£o||Ad||i,(^). 

Choosing £q > 0, such that Sq < I//4096, we have 

j/it) + j (/i|Vv|2 + (A + /x)|divv|2 + ^|Ad|2^ dx 

< l2Sue{cm + 12)||do||^.(^) (l|do||^.(n) + ^) := 9i- 

By fimi) . we obtain 

l|VM||^.(,.) <8 ( '^"fJ^°^ 9||VMo||^.(.)), 

where 

llV'dllL(^) := 5^ IIM-rfdli^(n)- 
Moreover, recalling |d| = 1, (12.160 and (12.90 . we see that 

Finally, utilizing (12. 9p . (I2.17p . Holder's and Poincare's inequalities, we get from the equation 
(Oj) that 

\\dtdj\\LA/:,(i^L\n)) <l|v • Vrfj-||i4/3(j,i2(n)) + ^IIArfj + | Vdpdj ||i4/3(j,i2(n)) 

<Ci(l])||Vv||L2(Q^)||Vd,||L4(Q,)+eTV4||Ad,. + |Vd|2d,.|U2(Q^) 
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In addition, we can also deduce that 

\\dtdj\\L^i,{min))*) <ci(f2) 

where {H^{fl))* denotes the dual space of H^{fl). 

Summing up the above estimates, we conclude that 

Proposition 2.1. Let Q be a bounded domain. We have the following a priori estimate for the 
initial- boundary value problem with initial data |do| = 1.' 

SUp(||ypv||L2(n) + ||p||L7(f1) + ||Vd||L2(f^)) + ||Vv||i2(j^L2(j^)) < C2{So). (2 18) 

t&I \ ■ J 

Moreover, if the initial energy satisfies Sq < Z//4096, then 

||V^d||L2(Qj,) + ||Vd||L4(Qy) + ||f^td||i4/3(7^L2(Q)) + 1 1 C^t d 1 1 ^2 (/^(Hl (Q))* ) 

< C3(||do||H2(f7),i'o)- 

Here C2 and C3 are positive constants which are, in particular, nodecreasing in their variables. 
Moreover, 02 also depends on the given physical parameters A, fx, v and 9; and C3 on v, 9 and 
the domain VL. 




3. Strong solvability of sub-systems in the third approximate problem 

In order to get a weak solution to the problem fll.ip - fll.Sp in a bounded domain f2, we shall first 
investigate the existence of solutions to the third approximate problem of the original problem 
(0)-([I3]): 



dtp + div(pv) = eAp, 

a^d + v • Vd = ^(Ad + I Vdpd), 

{p^){t) ■ *dx - uiQ- *dx 



pAv + (p + A) Vdivv - AVp'^ - 5V p'^ - e(Vp ■ Vv 

VdPl 



Jn 

- div(pv (g) v) - z/div ^Vd ® Vd 

for all t G / and any ^ G X„, with boundary conditions 

Vp • n|an = 0, \\d^ = 0, d\QQ -- 
and modified initial data 

p(x,0) = po e 0<p<po<p<oo 
d(x, 0) = do G H3(fi), v(x, 0) = vo G X,, 



*dxds 



do, 



(3.1) 
(3.2) 



(3.3) 



(3.4) 

(3.5) 
(3.6) 



where n denotes the outward normal to dfl, and e,6,(3,p,p > are constants. Here we briefly 
introduce the finite dimensional space X^. We know from 3^, Section 7.4.3] that there exist 
countable sets 

{Xi}Zi, < Ai < A2 < ■ • ■ , and 

{^i}r=, C Wl'^iQ) n W^'^iQ), 1 < p < 00, 
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such that 

- /iA*i + A)Vdiv*i = Xi^i, i = 1, 2, ■ • • , 

and is an orthonormal basis in L^(f2) and an orthogonal basis in HQ(f2) with respect to 

the scalar product J^[fidjU ■ dj\ + (/x + A)divudivv]dx. We define a n-dimensional Euclidean 
space X„ with scalar product < •, ■ > by 

Xj = spanl^jj^'^i, < u, V >= / u ■ vdx, u, v G X„, 

and denote by ^„ the orthogonal projection of L^(i7) onto X„. 

In the next section, we will show the unique solvability of the third approximate problem based 
on the idea from [sj Section 7.7] with certain restrictions imposed on the initial approximate 
energy. Moreover, the unique solution enjoys the energy estimates such as (12.181) . (12.191) . As a 
preliminary, this section is devoted to the global solvability of the Neumann problem (13. ip for the 
density and the local solvability of the non-homogeneous Derichlet problem (13. 2p when v is given. 
Due to the difficulty of the supercritical nonlinearity |Vdpd, we need here higher regularity im- 
posed on the initial data do (see (13. 6p ) in order to get the local strong solution d G L'^il^y H^(fi)) 
for some time interval I^. Thus, using the embedding theorem, Vd G L^(J^, L°°(n)). Such 
regularity is very important to deduce the pointwise property of d (see (12. ip ). 

3.1. The Neumann problem for the density 
We consider the following problem: 

dtp + div(pv) = eAp, (3.7) 

with initial and boundary conditions: 

p(x,0) = po(x) G W^^•°°(^]), < p < po(a;) < p < cx), (3.8) 
Vp -1119^ = 0, (3.9) 

which can be uniquely solved in terms of v G L°°(/, WQ'°°(fi)). The global well-posedness of the 
above problem can be read as follows (see, e.g., |34, Proposition 7.39]). 

Proposition 3.1. Let < a < 1, Q be a bounded domain of class C^'", and po satisfy Ii3.8\) . 
Then there exists a unique mapping 

^,„:L-(J,W^'°°(fi))^ C\I,H\Q)), 

such that 

(1) S^pQ ( v) belongs to the function class 

Ut := {pIp G L^{I,W^^\^))nC^{I,W^^''{^)),dtp G ^2(1, L^(^])) } , 1< g < oo. (3.10) 

(2) The function p = satisfies (X?) a.e. in Qt, ^3. 8^ a.e. in Vt and ^3. 9^ in the sense 
of traces a.e. in I . 

(3) ^poiy) is pointwise bounded, i.e., 



pe 



-/o l|v|lwi.-(n)dr ^ ^^^(v)(x,t) < pe^o l|vWllwi.-(a)d-^ ^ ^ j ^^g^ ^ G fi. (3.11) 
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(4) -(f ||v||l-(/,wi>°°{q)) < i^v, then 

\\ypoM\\L-iiuHHn)) < Ci||pobHn)e^(^"+'^"^ (3.12) 

I|V2^p,(v)||l.(q,) < ^v^||pobi(n)«:.e^(-+''?)* (3.13) 

l|5t^Po(v)||L2(QO < CiVt||pobi(f7)«:.e^(«"+''')* (3.14) 

for any t & I, where It '■= {0,t), and Qt := Q x Jj. The constant Ci in \3. 12^ - ^3.14^ 
depends at most on Vt, and is independent of k^, e, T, po o-nd v. 

(5) ypoiy) depends continuously on v, i.e., 

ll[^po(Vl) - ypoiy2)]it)\\mn) < C2iK.v, £,T)t\\po\\Hi{Q)\\Vi - V2||Loo(j^_wl.°°(f2)), (3.15) 

ll^t[^Po(vi) - ^po(v2)](t)||L2(Qt) < C2iK,v,£,T)Vi\\po\\m(n)\\vi - V2||L-(/t,wi>°°(n)) (3-16) 

for any t G /, and for any \\vi\\L°°{i,w'i-^°°{n)) < and \\v2\\L°°{i,w'^^°°{n)) < i^v The constant 
C2 is nondecreasing in the first variable and may depend on Q. 



Proof. All the results above have been shown in 3J, Proposition 7.39], except fl3.16p . Hence it 
suffices to verify f l3.16p . For simplicity, we always use the positive constant C2(k„, e, T) to denote 
various constants depending on k^,, e and T. 

Let pi,p2 G TIt be two solutions of the problem f l3.7p -f lXI?]l with v replaced by vi and V2 
respectively. Then after a straightforward calculation, we find that 

\jt I |V(Pi-P2)rdx + e I |A(pi-p2)|Mx 

= y (^V2 • V(pi - P2) + (Pi - P2)divv2 + pidiv(vi - V2) (3-17) 
+ (vi - V2) ■ Vpi) (Api - Ap2)dx. 



Using Cauchy-Schwarz's and Poincare's inequalities, the right hand side of f l3.17p can be bounded 
from above by 

Ci(r2) /.. WTif \ll2 II ||2 II ||2 \ 
(^l|V2||wi.°=(n)l|V(P2 - Pl)|lL2(n) + 11^1 - V2||wi>°°(n)IIPlll//i(!^) ) 

+ ^||A(pi — P2)|Il2(j^) for some constant Ci{^l) depending on f2, 

which, together with (13.170 . yields 

d 
di 



< 



|V(pi-p2)rdx + £ / |A(pi-p2)rdx 



(l|v2||wi.°=(Q)ll^(P2 - Pi)\\h{n) + 11^1 - V2||wi.°°(n)IIPl|lffl(!^)) 



e 

Therefore, by Gronwall's inequality and (I3.12p . 

||V(pi - P2)(t)||L2(n) < C2{k^, e,T)Vi\\po\\m{n)\\^i - ^2\\L^{it,wi---{n))- (3-18) 
Moreover, we have 

||A(pi - P2)||l2(qo < C'2(/«^,£:,T)v^||po||//i(n)||vi - V2\\L°-{it,w^-°°(n))- (3.19) 
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Consequently, using fl3.18p . f l3.19p . we get from fl3.7p that 

\\dt{,pl- p2)\\LHQt) 

= II (v2 • V(pi - P2) + (Pi - P2)divv2 + pidiv(vi - Vg) 
+ (vi - V2) ■ Vpi) - sA{pi - P2)||l2(qo 

< Ci{Q) (^/t^||V(pi - P2)||L2(Qt) + V^||Pl||L°°(7t,Hi(f7))||vi - V2||L°°(7t,Wi>°°(n))) 
+ £||A(pi -P2)||l2(q,) 

< 6*2(^1,, £, T)\/t||po||jfi(n) II vi — V2\\L^{it,w^-°°{n)), 
which imphes (13. 16^ . 

3.2. The local solvability of the direction vector 

Now we turn to show the local solvability of the direction vector. More precisely, we will show 
that for 

V G V := {v I V G C°(/, U\n) n L°°(r])), dtv G l2(/, h^(i^))}, 
there exists a G (0,T], such that the equation 

9td + v Vd = ^(Ad+|Vd|2d) (3.20) 

has a unique strong solution d(x, t) on [0,T^) satisfying the initial and boundary conditions 

d(x,0) =do(x) G H3(^]), (3.21) 
d|an = do(x) fortG(0,T;). (3.22) 

This local existence can be shown by modifying the arguments in |J, Section 3]. Here we give 
its proof, since some arguments in our proof are very different from those in jl. Section 3], and 
can be applied to the Cauchy problem in Section [6] by using the Ehrling-Nirenberg-Gagliardo 
interpolation inequality and domain expansion technique. Moreover, we shall show that the 
solution d continuously depends on v, which will be also used in the proof of the local existence 
of solutions to the third approximate problem. 

3.2.1. Linearized problems 
Denote 

w := {b G c\i:,ii\n)) I d,h G c'{i:,L'{n))nL'{i:,Hl{n))}, 

V^:={vGV|||v||v<ir}, 

:= {b G D* I b(x,0) = do(x), ||b||p, < 

where K and are positive constants. 



I l|2 I II ||2 I \\f) ||2 

\^\\co{i,m{n)) + \\^\\co{i,L°°{n)) + \\'^t^\\L'2{i,H^{n)) 

l|t>||lD)* := ^l|f^tb||^2(/._Hi{f^)) + ll^tb||^cx)(7*^L2(Q)) + ||b|||oo(/*^H2(f7)) + l|b||L2(/*^H3(0)) 

and := (0, TJ) C /. Without loss of generality, we assume > 1. 



(3.23) 
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To show the local existence of solutions, we will construct a sequence of approximate solutions 
to fl3.20l) - fl3.22p and use the technique of iteration based on the set D*^. First we linearize the 
original system f l3.20p for b G D*^ as follows: 

dtd-eAd = e\Vh\%-vVh (3.24) 

with a given function v G Yk, and initial and boundary conditions 

d\t=o = do, X e n, (3.25) 
d\dn = do, t e II (3.26) 

Since F(b) := [Vb^b - v ■ Vb G L2(/*, H^^^)) and dtF{h) G L^I*, (U^n))*), we can apply 
a standard method, such as a semidiscrete Galerkin method in [6], to prove that the initial- 
boundary value problem fl3.24p - fl3.26l) has a unique solution 

d G cv:, ii'm n L'iii u'm, did g i^iii iu\n)n 

d,dec'{r,,VmnL'{r„iilm. 

Here we have used (H^ ((]))* to denote the dual space of H^(fi). Therefore, we can get a solution 
d^ to fl3.24p with b replaced by some given d'^ G D*^. Assuming d'^^^ G for > 1, we can 
construct an approximate solution d*^ satisfying 

dtd'' - 9Ad^ = 0|Vd^-^|M'^-^ - V ■ Vd'^-^ (3.27) 

with initial and boundary conditions: 

d^|t=o = do, x G 
d'^lan = do, t G II 

3.2.2. Uniform estimates 

Next we derive the uniform estimates: 

\\d%w < (3.28) 

for some constant depending on do and the given constant K. We mention that in the following 
estimates the letter c > 1 will denote various positive constants independent of k^, v, do and 
VL] the letter Ci(f2) will denote various positive constants depending on Q, and the letter C(. . .) 
various positive constants depending on its variables (it may depend on Vt sometimes, however 
we omit them for simplicity), and is nondecreasing in its variables. Of course, they may also 
depend on the fixed value 6. 

First, we can deduce from fl3.27p and integration by parts that 

~ I |d¥dx + ^ I iVd^dx 
2 dt 

= 9 [ iVd'^-ipd'^-i ■ d'^dx- [ V ■ Vd'^-i • dMx + [ n-Vdo-dods 
Jn Jq Jdn 



and 



1 d 
2~dt 



Vd^pdx + ^ / lAd^dx 

Jvt 

|y^fc-i|2^jfc-i . ^dMx + / V ■ Vd^-1 • AdMx, 

n 
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which, together with Cauchy-Schwarz's and Holder's inequahties, imply that 
1 d 



fe-ll ||2 



2^l|d1l^i(n) < (^^ll|Vd'=-^|||t.(^)|||d'=-^|||i.(o) + ll|v|||i.(^)|||Vd 
+ l|diL2(Q) +^l|Vdo||L2(9n)||doi|L2(an). 

Hence, using the embeddings H'^iVt) ^ for d, H^iVt) )■ L^iVt) for (Vd,v), and trace 

theorem, we have 

|i|d1|^.(^) < C{K, ||do||H^(n))(«:^ + kIK') + 2||d^||^,(^), (3.29) 

which, together with Gronwall's inequality, implies 

|ijjfc/+Mi2 ^ ii2 , r^/T^ \\A II _ \t.^& , ,.2r^2^^^„2^ 



Hi(n) < (l|do|lHi(n) + C{K, ||do||H^(f,))(/s:;^ + KlK'')t^ 
for any t e I^. Taking = k^^ < 1 and letting TJ < Tj", we conclude 

I|d1|^i(n) < l|do||H^(Q)). (3.30) 

Next we derive bounds on dtd. It follows from (13.271) and integration by parts that 

-— I |9fd¥dx + ^ ! iVSid'^lMx 
2 at 

= -2e( I (Ad^-i ■ dtd^-^)d^-^ ■ dtd^fb^ + I Wd\-^ ■ S7d]-^dtd^^-^dtd)(bL 

+ / dtd'^-'d';-^Vdt^ -Vdtd^d^ +9 [ iVd'^^^fdtd''-' ■dtd'^dx 
Jn J Jn 

V ■ VStd'^-i ■ (9idMx - [ (9tv • Vd'^-i ■ 9tdMx. 

Jn 

Noting that {dtd, dtv, v) = on dVL, applying Cauchy-Schwarz's and Holder's inequalities, and 
fl2.12p . the right hand side of (13.311) can be bounded from above by 



c((||Ad'=-l^.(^Jd'=-l^.(^) + ||Vd^-l^.(^))||9,d^-iL.(^)||V5,d^-i£.(^ 

+ l|d'^-1^o.(f,)||Vd'=-l2,(^)||9,d'=-lL.(n)||V9,d'=-lL.(n) 



+ l|v||^o.(^)||V5,d^-iL^(n) + ||Vd'=-1^4(^) 

+ (||Va,d^-iL.(o) + ||5tv||L.(n)||Va,v||L.(^))||5,d^||L(^)) + ^-\\Vdtd%l,^^y 

Using H'^iVt) L'^iVt) and H^iVt) L'^iVt) for d and Vd respectively, the above term can be 
further bounded from above by 

Ci(fi) {{k^ + K'){\\VdA'-'UH^) + II V9,d^-i|,(^)) + 

(3.32) 

+ (||V5,d'=-iL^(f,) + ||9iv||L.(o)||V9iv||L.(o))||atd'=||^.(^)) + ^-WVdA'Wh^ny 
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Hence, we have 
d 



Jd,d'^\\lHn)+0\md'^\\Mn) 



< an) (^{k' + K'){\\Vdtd'-'h.^n) + \\Vdtd'-Tl^n)) + (3-33) 
Consequently, using Gronwall's and Holder's inequalities, we get 

for any t G Now, taking = k^^^ < 1, letting < T2 , and noting that 

||5td(0)||L2(n) = l|^(Ado + IVdolMo) - Vo ■ Vdol|L2(n), 

we conclude 

sup ||9td''||L2(j^) < C{K, ||do||H2(n))- (3.34) 

Moreover, using (13.271) . (I3.33p . (I3.34p . Cauchy's inequality, and the elliptic theory, we obtain the 
following estimates, 

||5id'=|U2(,.,HMf^) < C{K, ||do||H2(n)), (3.35) 

72 Jfc||2 IIV72J ||2 , II Q jfc||2 



+ l|v-Vd'=-||^.(,) + |||Vd'=-Td'=-^||^.(,)) (3-36) 
<C{K, ||do||H2(n)) + c|| V Vd'=-i^,(^) + c|| |Vd'=-i|2d'=-12^(^), 

and 

73jA:||2 ^nlJ^ \\A II ^ llwSj II _ \ , „IIV7/,, V7jfc-lM|2 

fe-l|2jA:-lM|2 (3.37) 
!(/,*,L2(n)): 

A:||2 



l|V^d'=||i2(,.,L2(f,) <C(ir, ||do||H2(n), ||V^do||L2(n)) + C||V(V . Vd'^-^)||i2(,.,L2(o)) 

+ c||V(|Vd'=-ird'=-i)||i2(,.,L2(^)), 



where 

l|V^d''||^2(j. L2(Q)) := ^ ll^j^i^mC?/ llL2(/.^i2(Q)). 

l<i,j,«,m<2 



To complete the derivation of (13.281) . it suffices to deduce the uniform estimates of the last two 
terms in (13.361) and (13.371) as follows. 

Since d G ©*, in view of [3J, Lemmas 1.65, 1.66], we have 

d(t)-do= \ (9td(r)dr for any t e 

which implies that 

||d(t) - dollHi(n) < v^||5td|U2(,*,Hi(n)). (3.38) 

Using (I3.38p . we obtain 

l|v ■ Vd'=-1^.(^) <c||v||^^(^) (||Vd^-i - Vd^i^^(^) + llVd^i^^(^) 
<ci^'(i'^d+l|do|lHi(n))- 
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Recalling that < t < < k^^^, we have 

l|v ■ Vd^-i^.(^) < C{K, ||do||Hi(n)). (3.39) 

To estimates the last term in fl3.36p and the the last two terms in f l3.37p . we shall introduce 
the following well-known Ehrling-Nirenberg-Gagliardo interpolation inequality: 

Lemma 3.1. Let Q be a domain in satisfying the cone condition. Ifmp > N, let p < q < oo; 

if mp = N, let p < q < oo; if mp < N, let p < q < p* = Np/{N — mp). Then there exists 
a constant c\ depending on m, N, p, q and the dimension of the cone C, such that for all 

||w||L9(n) < ci||t;||^^m,p||t'||^;"f2)' (3-40) 

where a = (N/mp) — (N/mq). 

Proof. The proof can be found in [H, Chapter 5]. □ 

In particular, we can deduce from fl3.40p that 

II^IIl-cc) < Ci{n)\\v\\H^(^n)\\v\\L2(n) (3.41) 
\\v\\%^^^<c,{Q)\\v\\H^n)\\v\\L^in) (3.42) 

for some constants Ci{fl) depending on Q. 

Making use of f l3.38p . fl3.4ip and f l3.42p . it is easy to see that 



<c^||Vd'' ^11^4(^)11^'' ^||L°°(n) + II ^11^4(^)1 



ifc-l||2 



\\L°°{n) 

llL4(f7)ll'-iO llL°°(f7) 



+ \\Vd'^-'-Vdt'\\iHnMr'"' 
<C^{n)(^K^d'-' - dollHMo) + l|do||^2(f,)^ 

for any t E I^. Since < k^^^, we have 

|||Vd^-i|M^-iL.(f,) < cm (l + ||do||^.(^)) . (3.43) 

Using Holder's inequality, f l3.4ip and f l3.42p . we infer that 
||V(v ■ Vd'=-^)||i.(,.,i^.(^)) + c||V(|Vd'=-i|M'=-i)||i.(,,,i^.(^)) 

< c^||Vv||^oo(j*^L2(f^))||Vd'' ^||i2(j*^Loo(j^)) + ||v||Loo(Qy)||V^d'' ^\\L2{i*,L'^{n) 

+ l|d^-i^o.(Q,)|||Vd'=-i|VM'=-li.(,.,L^(^)) + / ||Vd^-i^.(^)||Vd'=-l^.(^)ds 

< Ci{Q)(^K^\\Vd'''^\\Lo.(^i*,LHn))\\d'''^\\L^i*,HHn)) + K'^\\'^d'''^\\L^{i*,mn))Td 

II J A;— 1 II 5 II J k—1 1 II 

+ 11^ llL°=(/*,H2(n)ll^ \\\L^i2,n^n)) 
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which, together with the fact Tj' < k^^^ < 1, yields 

II V(v • Vd'^-i)||i.(,.,L^(^)) + c||V(|Vd'^-^pd^-^)||i.(,.,L^(^)) < + K^). (3.44) 

Finally, putting f l330|) . (ISlMl) - (13371 ). (EUS]) and (l3:ii|) together, we then obtain that 

lld'^lln. <C(ir,||do||H2(f,),||VX||L2(n)). 
Choosing Kd > max{C{K, ||do||H3(n), || VMo||L2(n)), ||do||H3(n), 1}, we get fl3.28p. 
3.2.3. Taking limits 

In order to take limits in fl3.27p . we shall further show that {d'^}^]^ is a Cauchy sequence. 
To this end, we define 

d^+^ = d'^+i - d'^ 

which satisfies 

dtd''+^ - Md'^+i = e[Vd^ : (Vd'^ + Vd'=-i)]d'= + e\V d''-^\'^d^ - v ■ Vd*^ (3.45) 
with initial and boundary conditions: 

d^+i|i=o = 0, xG^], 
d'=+^|an = 0, tell. 

Using Cauchy-Schwarz's and Holder's inequalities, (13.411) and (13.421) again, it follows from 
(I3.45P and the boundary condition that 

^P'^^i^Mo) + 2 (l|Vd'^+iL(c) + l|Ad^+1L(^)) 
<c(||Vd'=||^.(^)||V(d'= + d'^-^)||ioo(^)||d'^||i.(^) 

+ l|d^llL4(n)l|Vd'' ^||L°°(n)ll Vd''||L4(!^) + ||v^||Loo(Q)||Vd''||L2(Q) j + Hd'^+^l 1^1(0) 
<Ci(fi)||d'=||^.(^) {k^ + «:^(||d^||H3(^) + ||d^-iH3(n))) + l|d^+l^i(n)- + lld'^-lnscn)) 
Then, using Gronwall's and Holder inequalities, we have 



d 
dt 



d'^'\\l,r^^n)<Cim{Kh + Kl A||d'=||H3(0) + ||d'=-iH3(n))ds') sup ||d'=||^i(^) 
V ^0 / 0<s<t 



<C{K,Kd){t+Vi)e' sup ||di^i(^). 

0<s<t 



Taking 

T; = mm{T;, C-'iK, Kd)/iQAe'), 1/2} 
and letting ^T^, we have 



sup||d'^+i^.(^) < ^suplld'^ll^x(n). (3.47) 



Furthermore, from (I3.46p . (I3.47P and the elliptic estimates we get 



|d''^^||L2(7.^H2(Q)) ^ ^ sup ||d^||Hi(j^). 
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By iteration, we have 

sup lld'^+inMH) + l|d''+ii2(j*,H2(n)) < ^^^P PIIhmq) any k>2. 
In particular, 

oo 

(P1|L-(7-,Hi(n)) + ||d''||L2(^*^H2(n))) < OO, 

k=2 

which imphes {dfc}^^ is a Cauchy sequence in L°°{I}, H^(fi)) fl H^(fi)). Therefore, 

^ d strongly in H^(l])) n L\I*, H\Q)). (3.48) 

As a consequence, d is accurately a solution to problem fl3.20p - fl3.22p satisfying the regularity 
||d||D < Hd, by virtue of f l3.28p and fl3.48p . the lower semi-continuity, the elliptic estimate and the 
compactness theorem with time (see 0, Lemma 2.5]). 

3.2.4- Continuous dependence 

Finally we show that the solution d continuously depends on v. Let G (0, T|], vi, V2 G V^, 
di and d2 be two solutions of the initial-boundary value problem fl3.20p - fl3.22p . corresponding to 
V = Vi and v = V2 respectively. Moreover, two solutions satisfy 

||di||D* + ||d2||D* < oo. 

Multiplying the difference of the equations for di — d2 by di — d2 and integrating by parts, 
we obtain 

/ |di-d2|Mx + 2^ / |V(di-d2)|Mx 

= 2^ IVdiHdl - d2|'dx + j (V(di - d2) : V(di + d2))d2 ■ (di - d2)dx^ 

- 2 / (vi - V2) ■ Vdi ■ (di - d2)dx - 2 / V2 ■ V(di - d2) ■ (di - d2)dx. 

Similarly to that in the derivation of fl3.46p . we use f l3.4ip and fl3.42p to see that the right-hand 
side of the above identity can be bounded from above by 



C(||di||iD,*, ||d2||lD..,ir) (1 + ||di||H3(f7) + ||d2||H3{n)) 

X (j\di - d2||L2(Q) + ||vi - V2||L2{n)||di - d2||L2(n)) + 6'||V(di - d 

Hence, 

^l|di - d2||L(^,) + e||V(di - d2)||^2(f,) 
< C(||di||e*, ||d2||D*,^) (1 + ||di||H3(Q) + ||d2||H3(n)) 

X (^\\di - d2||L2(f7) + ||vi - V2||L2(f7)||di - d2||L2(n) 

in particular, 

— ||di - d2||L2{Q) <C(||di||in,*, ||d2||D*, -f^) (l + ||di||H3(f7) + ||d2||H3(n)) 

X (||di - d2||L2(f7) + ||Vi - V2||L2(n)) • 
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2;ilL2(n)' 



Thus, applying Gronwall's and Cauchy-Schwarz's inequalities, we conclude 

||(di-d2)(t)||L2(f,) < (t+ V^)C(||di||D.,||d2||B.,i^,T;)||vi-V2|Uo.(,^,L2(f,)) (3.49) 

for any t E I*^. Moreover, 

||V(di - d2)||L2(Q0 < (t + V^)C(||di||D*, ||d2||D*,^,T;)||vi - V2||L^(/,,L2{n))- 

Obviously, the uniqueness of local solutions in the function class D* follows from ( ]3.49p im- 
mediately. In particular, ||di||D*, ||d2||D* < i^d- As the end of this subsection, we summarize our 
previous results on the local existence of d. 

Proposition 3.2. Let K > Q < a < 1, VL he a hounded domain of class C^'", v G V^, and 
do G H'^(r2). Then there exist a finite time 

:= hi{K, ||do||H2(n), ||VMo||L2(f7)) e (0,min{l,r}), 

and a corresponding unique mapping 

^do ■ ^xivJith m place ofT) , H2(n)), 

where hi is nonincreasing in its first two variahles and Qt^ := Q x := Q x [0,T^), such that 

(1) ( v) helongs to the following function class 

n^K := {d I d G L\I^, did G L\I^, 

9,dGC°(/^L2(^]))^L2(/,^,H^(^]))}. 

(2) d = ^^,(v) satisfies BJW a.e. m Q^k , [IM) m Q and [IM) m Jf . Moreover, d\t=o = 
do. 

(3) S'^^iy) enjoys the following estimate: 

|K(v)||d* < C{K, ||do||H2(Q), ||VX||L2(n)), (3.51) 

where C is nondecreasing in its variahles, and in the definition ofW should he replaced 
hy . Moreover ( in view of h2.1\) and l{2.3\) ), 

\d{t) \ = 1 for any t, if |do| = 1; 

d2it) > do2 (^^y ^1^02 > ^02- 

(4) ^doi^) continuously depends on v in the following sense: 

ll«(vi) - K>M]{t)hHn) + ||V«(vi) - ^£(v2))||l.(qo 

r 3 (3.52) 

< VtC{K, ||do||H2(n), ||V do||L2(f7))||vi - V2||l<^(/j,l2{q)) 

for any t G , where C is nondecreasing in its variahles. 

Remark 3.1. We can further show that ^^,(v) continuously depends on v as the form of f l3.16p . 
However, the estimate fl3.52p is sufficient to prove the local existence of solutions to the third 
approximate problem. It should be noted that the constants hi and C above depend on the 
domain VL. However, in Section |6l we will see that the above result in bounded domains can be 
generalized to the Cauchy problem. 
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4. Unique solvability of the third approximate problem 

In this section, we establish the global existence of a unique solution to the third approximate 
problem f l3.ip -( l3l6l) . We first use the iteration technique and fixed point theorem to show the 
local existence. For this purpose, with the help of Proposition 13.11 and 13. 2[ we shall introduce the 
operator form of the approximate momentum equations (13.31) to construct a contractive mapping. 

4-1- Operator form 
Given 

peC\l,L\Q)), dtpeL^Qr), ess inf p(x, t) > > 0, ^i) 
we define, for alH G / 

by 

< ^p(t)V, w >= / p(t)v • wdx, v,wGX„. 
Jn 

Recall that any norms of X„ are equivalent, in particular, 

W^''P'{n) and (Wo^'^2(l]))*-norms are equivalent on X„, (4.2) 

where [Wq^'^^Q))* denotes the dual space of (Wo''P'(fi)), ki and k2 are integers, and < ^2 < 
oo, 1 < P2 < C)0, < A;i < 1 and 1 < pi < C)0 (or /ci = 2, 1 < pi < oo). Note that this property 
will be repeatedly used in the estimates below. First one has 

IIAwll^{x„,x„) < c(n) / p(t)dx, tel. (4.3) 

Jq 

It is easy to observe that -^pil) exists for all t G / and 

IK,(,')||if{x„,x„)<^, (4.4) 

where =Sf (X„, X„) denotes the set of all continuous linear operators mapping X„ to X„. 
By virtue of (14. 3p and ( 14. 4p . we have 

C\Tl) — 

ll^p7t)AiW^pwll-^(X".x.) < ^IIPiWIIli(q), t G /. (4.5) 

For the difference, the following inequality 

W^Piit) - '^P2WIU(X„,X„) < c{n)\\{p2 - Pi){t)\\Li(n), t G /, (4.6) 
holds. Due to the identity -^p^(i-) — -^p^it) ~ *^p^(t)('^pi(*) ~ ^p2{t))'^p^{t)i ih.aX 

\\-^;'it) - ^;'i^tM^M < ^ii(P2 -pi)(t)iiLi(^.), t G / (4.7) 

for pi{t),p2{t) satisfying ( 14. ip . 

Next, we shaU look for T* C (0,Tf'] and 



V G A := {v G C(/:,X„) I dtv G L2(/;,X„)}, := (CT^) C (0,T,^^ 
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with ||v||c(j.,H2(n)) + l|5tv||L2(/.,Hi(n)) < K for some K, satisfying 



(pv)(t) • *dx - / mo ■ *dx 

t 

' /iAv + (/i + A) Vdivv - AVp^ - 5Vp^ - ^(Vp ■ Vv) (4i 



Jn 



div(pv (g) v) - z/div Vd (g) Vd - — 



■ *dxds 



for all t G [0,Tn] and any ^ G X„, where p{t) = [^po('^)]('^) solution of the problem (13. 7p - 

(K9h constructed in Proposition [3H d(t) = ^£(v)(t) is the solution of the problem (I3:20|) - (l3:22|) 
constructed in Proposition l3.2[ By the regularity of (p, d) in Propositions 13 . 13 . 21 and the operator 
^p{t), the equations f l4.8p can be rephrased as 

= ^[So(v)]W {i^^o + ^[^(^p„(v), V, <(v))]d.) (4.9) 

with mo = (pv)(0), where := ^„ is the orthogonal projection of L^(i7) to X„, and 

^(p, V, d) = pAv + (/i + A) Vdivv - AVp^ - 5Vp^ - e(Vp ■ Vv) 

VdPA (4.10) 



- div(pv (g) v) - z/div ( Vd (g) Vd 

Moreover, one has 



aiv(t) =^[^^^^(^)](^) ^a,[^^^(v)](t)^[,^|^^ I ,?^mo + [^^(^po(v),v,^do(v))](s)ds 
+ ^So(v)]w[^-^(^Po(v), V, <(v))](t). 

(4.11) 

Auxiliary estimates 

We shall derive some auxiliary estimates on (p = ^pg(v),v,d = ^<^(v„)) and (p^ = 
^P(,(vfc), Vfc, d = ^^(vfc)), = 1, 2, where v and v^ belong to the class A, and 

l|v||c;(/*,X„) + ||5iV||L2(/.^x„) < K, \\Vk\\c{T*,X„) + ||5tVfe||L2(7*,X„) < K, 

with being a positive constant. By the equivalence of norms in (14. 2p . we have 



We denote (^(v), ^(v)) = (^po(v), ^£(v)) for simphcity. 
From ( KlUh and (jO]) we get 

||^^(p, v,d)||x„ <c(n)(||v||x„ + ||p|U-(n)(||v||x„ + llvll^J 

+ IIPllI°°(n) + IIPllL°=(f7) + ll^d|lL2(n)) • 



(4.12) 
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From fIXTT]) . flX^ . (ESD and KT^ . it follows that 

||^^(^(v), v,^(v))(t)||x„ </i2(if,P, ||do||H^(f,),||VMo||L^(n), T,n), t E i:, (4.13) 

where the constant h2 is nondecreasing in its first three variables. 

Employing flCTj) and the formula F{zi)-F{z2) = J^' F'{s)ds with F{s) = and F{s) = s^, 
we obtain 

< ^(pi, vi, di) - o/r(p2, V2, da), $ > 

= / [/iA(vi - V2) + (p + A) Vdiv(vi - V2)] • $dx 
Jn 

+ / [(Pl - P2)u\u{ + P2K - ^^2)^ + P24K - M2)]^j'^*dx 

+ e [ p2[A(vi - V2) ■ $ + - ui)dj^']dx + e [ (pi - P2) (Avi ■ # + 9,uia,$0dx 
Jn Jn 

+ u [ did2{djd'l - djd'^)dj¥dx + u [ djd\{did\ - 9i4))9j$Mx 
Jn Jn 



''"(7As^-^ + 5/3/-i)ds 

'n L"'p2 




div$dx 



(4.14) 

for any $ G Wo'°°(fi). Making use of flSTTT]) . f l33T|) . (1421) . f liAij) . and the elementary properties 
of the projection := (see js^, Exercise 7.33]), we get 

||[^^(pi,Vi,di) - ^^(p2,V2,d2)](t)||x„ 

< /i2(ir,p, ||do||H2(n), ||VMo||L2(n),T,n){||(vi - V2)(t)||x„ (4.15) 

+ ||(P1 -p2)(t)||Li(n) + ||V(di - d2)(t)||L^(n)}, t e 

where is again nondecreasing in its first three variables. 

Thanks to 04.71) . fl3.15p . (13.111) . and Holder's inequality, we obtain 

ll-^I^VoiW ^ -^[^V2)]wll-^(^"'^") - ||Polki(n),P,^,^)i||vi - V2||co(/,,x„), (4.16) 

for any Vi, V2 G C°(/*,X„), and t G /*. Here the constant /is is nondecreasing in its first two 
variables, and is nonincreasing in its third variable. Since 



/ V^(^(v), V, ^(y))\{s)ds - ^[:;^(.)](,,) / '[^^(^(v), V, ^{y))\{s)ds 
Jo Jo 

^[>(v)](*i) r [«^(V)^(^(V), V, ^(V))] (S)d. 



tl 

+ [AyUit.)-Ayi^m2)) [^^(^(v),v,^(v))](.)d., 

we make use of (EH]), (S3]), (|42]) and (14A31) to find that 



t2 



^GUt) I [^-r(^(v),v,^(v))](.)dsGC°(/:,XO. 



And notice that 

< h,{K,p,\\d4nHn),P,\\^'dohHn),T,n)t, (4.17) 

C(/t,X„) 
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^yU A^^(^(v),v,d](s)ds 

•/ 



for any t & I*, where /i4 is nondecreasing in its first three variables and nonincreasing in its 
fourth variable. Similarly, ^^^^^^^^^{^mo) E C(/*,X„) and 

||^^|^)(^mo)||c(/,x„) < p-^e^i^mo||x., t E C (4.18) 
Finally, we denote 

u = ^yl^^^a^yn^yl^) (^^mo + £[^^(^(v), v, ^(v))](s)d.^ 
+ ^-;^)[^^(^(v),v,i^(v))]. 
Thus, we can use (13.111) . (14. 4p . (14. 5 p and (14.130 to evaluate (I4.19p and obtain that 

||u||x„ < HK^p, ||do||H2(n),P, ||VMo||L2(c),T,ri)(||at^(v)||Li(n)(l + ||^mo||xJ + l), 
which, together with (I3.14p and Holder's inequality, yields 

||u||L2(/t,x„) < h{K,p, \\po\\m{n), ||do||H2(n),P, || V^do||L2{c), 7", ^)v^(l + || <^mo||x„), (4.20) 
where is nondecreasing in its first four variables and nonincreasing in its fifth variable. 
4-3. Local existence 

Now we are in a position to prove existence of a local solution by applying a fixed point 
theorem. To this end, we assume that 

5max| l'^^°l'^" , ||v(0)||x„} < (4.21) 

and take 

^ ^r-ll II lU II T<^^ ^ ■ ln(5/4) 25 ^\ 

To :=To(p,||po||//i(f7),||do||H3(c),P,i^,T,n) = mm| — ,— ^,^,^^-^,T^ |, (4.22) 

so that To is nonincreasing in its first three variables and nodecreasing in the fourth variable. 
With this choice we have 

/i4To<^, p-ie^^°||^mo||x„< J, /^s V^(l + || ^mo||xJ < y ■ 
Therefore, by virtue of KT7\f . fHlSD . flCTjl and fOTD . the mapping 
r : A ^ A, 

r(w) :=^j>i^^j{^mo+ r[^^(^(w),w,^(w))](s)d4, (4.23) 

t/ 

maps 

BK,ro = {w G A I ||w||c(/,^,X„) + ||<9iW||i2(j^^,x„) < K} 

into itself for any < To < To, where we can take X„ = HQ(^7). 

In the next step, we prove that T is contractive. Keeping in mind that 

- = i^-^ - ^p-^)(wi) + ^-^\w^ - W2), 
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we get from fl4.23l) that 
T(wi) - T(W2) 

+ -^yU) f (^^(^(Wl), Wi, ^(Wi))) - ^^(^(W2), W2, ^(W2)))) {S)ds. 

Recalling that p < p(0) < p, one has 

||^mo||x„ < cHp||v(0)||x„ < c{n)pK. (4.25) 

We apply KT3\f . f l4J6|) and g^S]) to bound the first term in ^Ml, and use fl3TT]) . f l3J[5|) . f l332|) . 
f l4.4p . ( I4.15P and Holder's inequality to majorize the second term in (14 .24^ . and obtain 

Iir(w0-r(w2)||x„(t) 

^ (4.26) 

< hfiip, WPoWmin), i|do||H2(Q),P, ||V do\\i^2(^^), K,T,n)Vt\\w:L - W2||co(j,,x„), 

for t G [0,To], Wi,W2 G Bk^To-, where is nondecreasing in its first three variables and nonin- 
creasing in the fourth variable. 

Similar to (14. 24 p . we also have the identity 

dtTiwi) - dtT{w2) 

= i-^yU) - + ^ [=^^(^(wi), wi, ^(wi))](s)ds 

+ ^^;w2)(^9..^'{wi) - ^a„r(w2))^rU){'^™0 + / [<^^(^(wi), Wi, ^(Wi))](s)ds 

+ J^y\^^^J^a,.n^,){J^y\^^^ - ^^;^^))|^mo + ^ [^^(^(wi), wi, ^(wi))](s)ds 



+ ^^(w2)^9.-^(w2)-^r(w2){ ^ [,i^^(^(Wi), Wi, ^(Wi)) 



- ^^(J-(w2),W2,^(w2))](s)dsj + (^^;^^) -^_^,;^,))[^^(^(Wi),Wi,^(wO)] 
+ ^^(w,)[^^(^(wi), Wi, ^(Wi)) - ^^(^(W2), W2, ^(W2))]. 

Now, employing flSll]) . fl332D . (gSD, (|13D, dMD, fHA3|l . fl4J5|) . fICTD and fOSD to control the 
six terms on the right-hand side of the equality above, we deduce that 

||Wwi)-Ww2)||x„(t) 

< ^6(P, IIPollHi(n), ||do||H2(Q),P, ||VMo||L2(n),^,T,n)(^||at[^(wi) - =5^(w2)](t)||Li(Q) (^4,27) 
+ (||9t^(wi)(t)|Ui(f,) + 1)||W1 - W2||cO(J„X,0 + l|V(di - d2)(t)||L2(n) 



for a.e. t G [0,To] and for any wi,W2 G Bk^Tq: where is nondecreasing in its three variables 
and nonincreasing in the fourth variable again. 

Consequently, substituting (ISlD . (KT^ . and (13321) into KT!^ . we find that 

II Wwi) - Ww2)||l2(j^,x„) 

(4.28) 

< hfiip, IIPollHi(n), ||do||H2(f7),P, II V do||L2(n),^,T,n)Vt||wi - W2||co(j,,x„)- 
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Adding f05]) to fOg]) . we finally get 

lir(wi) - r(w2)||x„(t) + II Wwi) - dtriw2)\\LHiM 

< hfiip, ||po||Hi(f7), ||do||H2(n),P, ||VMo||L2(n),-f^,7',n)t 

(||wi - W2||c;0(j-,x„) + ||5t(wi - W2)|U2(j-,x„)) 



for any t G [0, Tq] 
If we take 



T: = min <! To, ^ 

2/i, 



then T maps Bk,t* C C(/t*,X„) into itself and is contractive. Therefore, it possesses in Bk,t* 
a unique fixed point v which satisfies (14.8p . Thus, we have a solution (p = ^(y), v, ^(v)) which 
is defined in Qt* and satisfies the initial-boundary value problem fl3.ip - fl3.6p for each given n. 
Moreover, we see that T* has the form 

< T* = h{p, ||po||//i(n)5 ||do||H2(f7),p, || V^do||L2(f7), ^, T, ra) < To, (4.29) 

where h is nonincreasing in its first three variables and nondecreasing in the fourth variable. This 
means that we can find a unique maximal solution (p„,v„,d„) defined in [0,T„) x Q for each 
given n, where Tn < T. 

4.4- Global existence 

In order to show the maximal time Tn = T for any n, it suffices to derive uniform bounds for 
and ^n(p„v„). However, we need to impose an additional (smallness) condition on 
the initial approximate energy to get the uniform bound of ||d„||£oo(-/^ H2(n))- We remark that T^ 
in fl4.29p depends on || V^do||L2(n)! since we have used the fact d'^"^ = do on dfl and the elliptic 
estimate to obtain fl3.37p . and thus we need an auxiliary term || V'^do||L2(n)- Such auxiliary term 
will not changes, if only the initial data of d(x, t) changes and the boundary value of d(x, t) does 
not change. This is why we do not estimate the uniform bound of || V^d„||^oo(/^ l2(q)). 

For simplicity of notations, we denote 

(p, V, d, ^m) := {pn, v„, d„, ^„(p„v„)). 

We mention that in the following estimates the letter G{. . .) will denote various positive constants 
depending on its variables. 

First, we derive energy estimates similar to Proposition 13. II Differentiating (14. 8 p with respect 
to t, integrating by parts, employing (13. ip and the regularity of (p, v), we obtain 



d 
dt 



where 



Es{t)+ [ (p|Vv|2 + (A + p)|divv|2 + e5/3p'3-2|y^|2^ jx<z/ /" (Vd)^Ad • vdx, 
Es{t) = I (^^1{P>0} + Qip) + p^P') 



Noting that d belongs to the function class (I3.50p with T„ in place of Tf, d satisfies (12. 7p . and 
one further has 

^Ssit) + f (p|Vv|2 + (A + p)|divv|2 + ue{\Ad + iVdpdH + 66(3p^-^\\/p\^) dx < 0, 
dt Jn^ 

(4.30) 
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where 

E,{t) := E,{p, m, d) := (^-^l^^^o} + + j^P^ + j 

with m = pv. Integrating f l4.30p over (0,t), we get 

/ / (Ai|Vv|2 + (A + p)|divv|2 + z/0(|Ad+ iVdlMp)) dxds 

< £^5(0) := £^5(po, mo, do) < ^5,0, 

where £5^0 is a given positive constant. In particular, 

II VPv|U-(/„,L2(n)) < 2^5,0, (4.31) 
||Vd|Uoo(,„,L2(f,)) < 24o/^^, (4.32) 
l|v||L2(j„^Hi(n)) < c{yL)S5fl/ii. (4.33) 

With the help of fl4.31l) - fl4.33l) . we can deduce more uniform bounds on (p, v). Using (13.111) 
and fl4.33p . thanks to the norm of equivalence on X„ (see (14.21) )). we find that 

Gi{p, Ss,o, T,n)<p< G2{p, Ss,o, T, n), (4.34) 

from which, f l4.3ip and f lO]) . it follows that 

I|v||c70(/„,x„) < G{p, Ss,o, T, n) (4.35) 

and 

ll<!^(pv)||co(/„,x„) < c(r2)||pv||co(4,L2(n)) < G{p,£s,o,T,n). (4.36) 
Applying fOSD to IKT2\i . one gets 

\\p\\co{i,„m(n)) < G{p, ||po||Hi(c),^<5,o,r, ""-)• (4.37) 
Utihzing f02|) . f011) - f06|) . arguing similarly to that for fOOj) . we obtain from fHTTD that 
||5tv||L2(/„^x„) < G{p,p, \\po\\m{n), \\do\\u^Q),^5fi,T,n), (4.38) 

Hence, we have shown the uniform boundedness of p, p, ||p||co(7„,x„)) ll'^(pv)||c'0(/„,x„)5 
l|v||co(J„,x„) and ||9tv||i2(/^ x„)- It remains to show the uniform boundedness of ||d||^oo(/^ H2{n))5 
which can be obtained by following the spirit of proof in Section 13.2.21 For the reader's conve- 
nience, we give the proof in the following. 

First, we rewrite the energy inequality fl4.30p in the following form as in (12. lip : 



4 

L4{n)- 



S{t)+ [ [ [p\Vv\^ + {X + p)\diYv\^ + iye\Ad\^ + e6(3p^'^\Vp\^]d^<S5,o + iye\\\Vd 
Jo Jn 

Let Ss,o = z^/4094. Arguing in the same manner as in the derivation of (I2.19p . we deduce that 

||VM||L2(Qy) + ||Vd||L4(Qy) + ||<9td||i4/3(/,L2{Q)) < ^ ( 1 1 do 1 1 H2 (C) , ^5,o) • (4.39) 
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Now we proceed to derive uniform bounds of higher derivatives on d. Differentiating 
with respect to t, multiplying the resulting equations by dtd in L'^{fl), recalling |d| = 1, we 
integrate by parts to infer that 

= 2 [ dtd - {OAdtd + ^1 Vdl^^td + 2^(Vd : Vdtd)d - S^v ■ Vd - v Vdfd) dx (4.40) 
Jn 

< -e\\vdtd\\l^^^ + (100 + i)|||9id||vd||i.(^) + 2||9*v||^,(^) + ^Ml^t^Qj\dtd\\l^^), 

where the second term on the right-hand side of f l4.40p can be bounded as follows, using f l2.12p . 
Holder's and the triangle inequalities. 

(io0 + i)|||a,d||vd||i.(^)<(io0 + i)ll|9,d|||i.(^)|||vd|||i.(^) 

2 

<{10e + l)Y,\\dtd4lHn)\\\^mlHn) 

i=l 

2 (4.41) 

<(io0 + i)^v^l|9tc?,|U2(o)||vaici,|U2(n)|||vd|||i4(o) 

i=l 

<^l|V9,d||^.(^) + ^^^^f^||9,d||^.(^)l||Vd|||t4(^). 



Inserting f l4.4ip into (14. 40 p . we conclude that 

|l|5*d||^.(^) + ^||V9,d||^.(^) 
(100 + 1)2 , 1 



< 



9 ^ 9 

Thus, applying Gronwall's inequality, we have 



ll|Vd|||i.(^) + -||v||i.(Q^) l|5td||^.(^) + 2||9,v||i.(^). 



^IIIVd|||!,_,+f||v|i^ 



l|5.d||i.(^)<(^||9,d(0)||l.(^) + 2||9,v||^.(Q^)je - — ".^(Q.r^-"™.). 
Noting that 

\\dtdiO)hHn) = l|0(Ado + |Vdo|Mo) - vq ■ Vd^h^^n), 
we use fOSj) . fOSj) and f lQQj) to arrive at 

+ liva,d|| 

Recalling that |d| = 1, it follows from (13. 2p that 

9^ [ |Ad|2dx<302 f |Vd|Mx + 3 [ |9td|Mx + 3 [ Iv^VdlMx. (4.43) 
Jn Jn Jn Jn 

Similarly to the derivation of fl2.15p . the first term on the right-hand side of f l4.43p can be 
estimated as follows. 

39^ I |Vd|Mx 

2 (4.44) 
< 38402(c(fi) + 12)ido||^.(^) (l|do||^.(o) + ^) + ^^^||Ad||i.(^). 
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Putting (gaSD and fICTD together, using (gSSD, (i32D and f lO^ . we get 
^ ^ I Adpdx <{c{n) + 12)ldo||^.(^) (l|do||^.(n) + ^ 



+ 3||c?td||L2(n) + 6||v||L^(f^)||Vd||L2(Q) 
<G{p, p, \\po\\m(n), \\d.o\\u2{n), £s,o,T,n). 



(4.45) 



Hence, by the elhptic estimate fl2Hl) . f 02|) . fH:i5|) and the fact |d| = 1, 

||d||Loo(/„,H2(n)) <G{p,p, \\po\\m(n), \\dQ\\u^n),£s,o,T,n). (4.46) 

The inequahties f l4.34p . f l4.35p . fl4.37p . f l4.38p and f l4.46p furnish the desired estimates which, 
in combination with f l4.2ip . f l4.22p and fl4.29p . give a possibihty to repeat the above fixed point 
argument to conclude that T„ = T, and the global solution (p„, v„, d„) is unique. To end this 
section, we summarize our previous results on the global existence of a unique solution (p^, v„, d^) 
to the third approximate problem f l3.ip - fl3.6p as follows. 

Proposition 4.1. Let 

5 > 0, P > 0, e > 0, and < p< p < oo. (4.47) 

Assume that fl is a bounded C"^'" -domain (a G (0, 1)), and the initial data (po,mo,do) satisfies 

Ssipo, mo, do) < £5,0 ■■= Z//4094, (4.48) 
0<p<Po<P, PoeW'^in), |do| = l, (4.49) 
vo G X„,, do G U\n). (4.50) 

Then there exists a unique triple {pn, v„, d„) with the following properties: 
(1) Regularity. 



Pn satisfies the regularity as in Proposition \3.1\ with T in place of , 

v„ G C°(/,X„), 9iV„ G L2(/,X„), Vp„ G L\l,El{n)), p„v„ G C\l,El{n)), 



(4.51) 



d„ satisfies the regularity as in Proposition \3.S\ with T in place ofTf. 

(2) (p„,v„,d„) solves h3.1\) and h3.2\) a.e. in Qt, and satisfies h3. 3|) and (p„, v„, d„)|t=o = 
(po, vo,do). 

(3) Finite and hounded energy inequalities: 

j^S^it) + F'^it) + ^£5/3p^2|Vp.|2(^)dx < m V'il), (4.52) 

and 

^n^) + ^*(-^"(^) + ^^'^/3p^'|Vpn|'(s)dx^ds<£5(po,mo,do) a.e. in I, (4.53) 
where J^'^it) := J'(p„, v„, d„) and S^{t) := £s{pn, m„, d„) with m„ = p„v„. 
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(4) Additional uniform estimates. 

|d„| = 1 m Qt, (4.54) 

||VM„||l2(Qj,) + ||Vd„||i4(Q^) + ||<9td„||i4/3(/x2(Q)) < Gdldollnacn)), (4.55) 

y^l|Vp„||L2(Q,) < 0(5), (4.56) 

\\Pn\\^^^^^^<Gis,6) (4.57) 



(^see Section 7.7.5.2] for the proof of li4-56{ ) and Iji4.5'7\ )), where G is a positive constant 



which is independent of n and nondecreasing in its arguments. Moreover, if e is not explicitly 
written in the argument of G, then G is independent of e as well. 

Remark 4.1. Here we have used the notation Eq{Q) concerning the spaces of vector fields with 



summable divergence introduced in 34, Section 3.2]. For the reader's convenience, we give the 
detailed definition. We set 

E5'P(n) = {g G iL''in)f I divg e LP{n)} ^ L'^in) 

with the norm 

l|g||-B9>P = llgllL'J + ||divg||LP(n)- 
For the sake of simplicity, i?P'P(f2) is denoted by Ep{Q). Then we define 

-El'P , „ ; -EP 



E^o'^{Q) = V{Q) and = V{n) 
The notation Eq^{Q) will be used in Proposition 15.11 

5. Proof of Theorem 11.11 

Once we have established Proposition 14. we can obtain Theorem ILl lbv the standard three- 
level approximation scheme and the method of weak convergence as in |8|, |3l| for the compressible 
Naiver-Stokes equations. These arguments have also been successfully used to establish the 
existence of weak solutions to other models from fiuid dynamics. Here we briefiy describe how 
to prove Theorem 11.11 and state the existence of solutions to the first and second approximate 
problems for the reader's convenience, we refer to 38| or jsl,!^^] for the details of the limit process. 



In view of the proof in [38|, it suffices to analyze the convergence of the supercritical non- 
linearity |Vdpd. In fact, using the uniform bounds in (I4.54P and ( 14.55P on d„, applying the 
Arzela-Ascoli theorem and Aubin-Lions lemma, and taking subsequences if necessary, we deduce 
that 

d„ ^ d strongly in C%I, L'^iQ)) n LP{I, H\Q)) n W^'P{Q)) (5.1) 
for any p > 1 and r G [1, 2), which implies that 

|Vd„|M„ ^ iVdpd weakly in L\Qt). (5.2) 
Moreover, using Vitali's convergence theorem, and taking subsequences if necessary, we obtain 
I Vd„pd„ ^ I Vdl^d strongly in UiQr) for any r e [1, 2). (5.3) 
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In addition, we also have the regularity d G L'^{I, H^(f2)) and dtd G L'^{I, (H^(f2))*). In view of 



34 Proposition 7.31], we get consequently 

d G C%I,H\Q)). 



(5.4) 



Putting (15. 2p and (15.41) into the proof of [381, Proposition 4.1], we immediately obtain a weak 
solution of the second approximate problem which is the weak limit (p£,V£,de) of the solution 
sequence (pn, d^,) as — )■ oo constructed in Proposition 14.11 Moreover, refereing to the 
conclusions in |3J, Proposition 7.31], the existence of solutions to the second approximate problem 
reads as follows. 

Proposition 5.1. Let 6 > 0, P > max{7,8}, e > 0, and < p < p < oo. Assume that Q is a 
bounded C'^'°' -domain (a G (0, 1)), and the initial data (po5mo,do) satisfies fi4-4S\ ), fi4-4^ o,'>T'd 



mo G Vin), do G B.\n). 

Then there exists a unique triple (p^, V£,de) with the following properties: 
(1) Regularity. 

>Oa.e. mQT, Pe e C%I , L^^^^m C\I , (Q)) , 1 < P < /S, 
pi G L\l,H\n)), dtPe G L-^{Qt), VV. e L^(fi), 



(5.5) 



10/3-6 5/3-3 



vp., p,v,gl2(/,Eo^^+^' m, 



Pedx 



podx. 



2/3 



tm,:=p,v, gC°(/,L^:4(1])), d,GC°(/,Hi(fi)), d^ = do G C°'i(a(]). 

(2) Weak-strong solutions. 

dtPe + div(p,v,) - eApe = m P'(Qt), (5.6) 
dt{pe^e) + 5i(p.v,t;^) - pAv, - (p + A)Vdivv, + VApJ + 5Vp^ 



+ z/div Vd. (g) Vd. 



iVdJ^I 



+ £Vp,-Vv, = zn(I)'(QT))' 



dtd, + V, ■ Vd, = ^(Ad, + iVd.pd,), a.e. m Qt- 



(5.7) 
(5.8) 



Moreover, (p^, v^, dj(x, 0) = (po,vo,do). 

(3) (pe,Ve,de) Satisfies the finite and hounded energy inequalities as in ^4.52^ and ^.5!^ , and 
the uniform estimates ( [^ . 54^ and ( [^.55| ). 

(4) Additional uniform estimates. 

WPe^eW 2^ + WPe^eW ^ + ||PeVe|| 103-6 

L°°(/,LWT(n)) L2(/,L^(n)) L^?+3^(Qt) 



6,9 + llpslli^i+i/n^) + ell Vpell 10/3-6 



+ llPelVe 
+ V^ll Vpe||L2(QT) + Vpe ■ Vv, 



£11 5/9-3 ^ C('^); 



where G is a positive constant independent of e. 
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Here we explain how the assumptions fl4.5Up in Proposition 14.11 become f lS.Sp in Proposition 
15.11 First, we have Proposition 15.11 for fixed initial data satisfying fl4.50l) . If they satisfy the 
assumption (15. 5p only, then we can approximate (po5nio,do) by a sequence (pom, mom, dom) 
which satisfies f l430|) . (|13HD, and 



EsipQm, mom, dom) Ssipo, mo, do) as m +oo; 



(5.9) 



moreover, dom — )■ do strongly in H^(i7)nC°'"'^(0). Then it is sufficient totake Vom = -Pn(mo/Po) ^ 
X„. Obviously, Proposition 15.11 holds with these new initial data, and those uniform estimates 
in Proposition 15.11 are independent of m. Therefore, repeating the limit process as m — )■ oo, we 
obtain Proposition 15.11 as well. 

Now, with Proposition 15 . 1 1 in hand, following the proof of 38|, Proposition 5.1] or SJ, Propo- 
sition 7.27], we can obtain a weak solution {ps,vs,ds) of the first approximate problem as the 
weak limit of the sequence (p£,V£, d^) as e — i- constructed in Proposition 15.11 Thus, we have 
the following existence result. 

Proposition 5.2. Let S > 0, f3 > max{7,8}. Assume that Q is a bounded C'^^'^-domain (a G 
(0, 1)), and (po,mo;do) satisfies M.9\) - n~Al\) with /3 in place of'-/, and l\4-4^ - Then there exists 
a unique triple {ps, v^, d^) with the following properties: 

(1) Regularity. 

PS e C%1 Ll.^^m n C°(/, L^m n L^+\R' x /), i < p < 

Ps > a.e. in Qt, ps = 0, = in (M^\f2) x /, 
< Psl^sl e L\l,Lvm(M?))nL\l,L7+-s(R^)), 



ms := ps^rs G L^(/,L 



6/3 
6+/3 I 



^))nc°(/,L:-(fi)), 



[dseC'iI,U\n)), ds\9n = do. 

(2) {ps,^s,ds) solves / l5.6]) -/ r575|) with e = 0, and also satisfies the finite and hounded energy 
inequalities as in ^4-5^ and ^4-5^ with e = 0. Moreover, (p5, v^, d5)|(=o = (Po; Vo,do). 

(3) For any b satisfying (TJ^ and [Tm) . the function h{ps) is m C°{I, {n))nC'^{I, LP{n)) , 
1 < P < . Moreover, 

dtb{ps) + div[6(p5)v] + [psb'ips) - b{ps)] divv = m V'{R^ x /). 
For any bk (k > 0) defined by 



bk{s) 



bis) ifsE[0,k), 
b{k) if s E [k, oo) 



b'{s) ^/se[o,fc), 

if s E [k, oo) 



{bu)'M 

with b satisfying (TTT^ and U.18\) . bk{ps) belongs to C°(/, L^(i7)), 1 < p < oo. Moreover, 
OMps) + div[6fe(p5)v] + [psibkY+ips] - bkips)] divv = m V'{R^ x /). 
(4) {ps,"vs,ds) satisfies the uniform estimates ( 4-54 ) and ( (^.55 ). Moreover, 
\\ps^&\ 
+ \\Ps\^5\ 
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L°°(/,LWT(n)) 
|2|| 



L2(j-^LWH(Q)) 



L'2(I,L^{n)) 
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for any constant 9 E (0,7 — 1) (the uniform estimate on \\p5\\L^+e(Qj.) can be found in 118] for the 
two-dimensional case). Here G is a positive constant independent of 6. 

Remark 5.1. We assume that the initial density po satisfies fl4.49p in Proposition 15.11 However 
the condition f l4.49p on po can be relaxed by the condition "po G L^(fi) and po > 0" in Proposition 
15.21 T his process can be dealt by the approximation method (similarly as in (15. 9p ). please refer 



to [3^, Section 7.10.7] for the detailed proof. 



Finally, with the help of Proposition 15. 2^ we can follow the proof of [38|, Theorem 2.1] or |3J, 
Theorem 7.7] to obtain a weak solution (p, v, d) of the original problem (ll.ip - (ll.4p which is the 
weak limit as 5 —> of the weak solution sequence (psyVsyds) constructed in Proposition 15.21 
This completes the proof of Theorem 11.11 

6. Global existence of large solutions to the Cauchy problem 

In this section, we will briefiy describe how to prove Theorem 11.21 on the Cauchy problem by 
modifying the proof of Theorem 1 1 . 1 1 and applying the domain expansion technique. 

6.1. Local solvability of the Cauchy problem on the direction vector 

First we establish the local existence of solutions to the following Cauchy problem on the 
direction vector: 

9td + v- Vd = ^(Ad+ |Vd|2(d + e2)) in x J (6.1) 

with initial data 

d(x, 0) := do e H2(M2)^ |do(x) + esl = 1, c?02(x) + 1 > ^02 > (6.2) 

for some given constant rfgg, where do2 denotes the second component of do, and the known 
function 

V e V := {v G C°(/, U\R^) n L°°(R2)) I dtv G L^I, ii\R'^))}. 

The local existence can be shown by following the proof of Proposition 13.21 and applying the 
domain extension technique. Next, we briefiy describe the proof. 

By virtue of (16. 2p . we have a sequence of functions {d°(x)}„>i C H^(]R^), such that 

|d°(x) + eal = 1, c/°2(x) + 1 > ^02/2, d° ^ do in H^R^) 

and suppd° C -B_r„ := {x G | |x| < _R„} for some Rn > 1. In view of Proposition 13. 2[ there 
exists a unique local strong solution d to the following Direchlet problem: 

atd + V- Vd = ^(Ad+ iVdpd) inBR^xI, (6.3) 

with initial and boundary conditions 

d(x,0) = d° G H^BrJ, d{^,t)U^^ = for any t G /. (6.4) 

Moreover d2 > 0(^02)- However, the constants hi and C in Proposition l3.2l depend on the bounded 
domain Q = -B_r„, since we have used the interpolation inequalities (I3.4ip and (I3.42p (of course, 
the embedding theorems H'^{n) ^ L°°{n) and H\n) ^ L^(fi) used in and (13:^ can be 
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replaced by such two interpolation inequalities, respectively), where the interpolation constants 
depend on the bounded domain Q = Br^. Fortunately, we have the following results: 



< cMU^iBnJHL^Ba^) and MW^.^b^j < c||t; j || ^^(b^ j (6.5) 

for some constant c independent of i?„ > 1, which can be deduced from f l3.4ip and f l3.42p with 
Bi in place of VL by scaling the spatial variables. In addition, the elliptic estimate used in fl3.37p 
is replaced by the following special eUiptic estimate: if G Hq{Br^), f G H^{Br^) and At; = / 
in V'{Br^), then 

) ^ '^\\f\\H^iBR„) for some constant c independent of i?„ > 1, (6.6) 

which can be deduced from the standard elliptic estimate on the domain Bi by scaling the spatial 
variables. With these facts, repeating the proof of Proposition 13. 21 with a trivial modification, we 
can obtain Proposition 13. 21 with Q = B^^. Moreover, the constants hi and C can be independent 
of Br^ and || V^do||L2(_Bfl )• Then, using the domain expansion technique (see jij), we can obtain 
a unique strong solution d of the Cauchy problem (16. ip . (16. 2 p as the weak limit of the sequence 
d_R„ with initial data d/j^|t=o = d° as Rn — )■ oo. Thus, we have 

Proposition 6.1. Proposition \3.2\ holds with in place of Q, where d = ^J^, (v) solves the 
Cauchy problem h6. I{6.^) . the constants hi and C only depend on K and ||do||H2{R2)) O'nd 

rf2(x,t) + 1 > ^02/2 for any (x,t) G x jf , «/c/o2(x) + 1 > ^02- (6-7) 

6.2. Global solvability of the approximate problem 

Once we have established Proposition 16.11 we can see, in view of the proof in Section HJ that 
there exists a unique local solution to the following approximate problem to the original problem 
([IIID-([I3D with VL = M^: 



/ (pv)(t) ■ *dx - / mo • *dx 

J Bn J Bn 



/iAv + (/i + A) Vdivv - AVp^ - 5Vp^ - e{Vp ■ Vv) - div(pv ® v) (6.8) 



'0 JBr 

- z/div ( Vd (g) Vd 



|Vd|2l\ 

*dxds for alH G / and any * G 



2 

dtp + div(pv) = eAp in Br x /, (6.9) 

^td + v- Vd = ^(Ad+|Vd|M) in x J, (6.10) 

with boundary conditions 

Vp ■ n\gBR = 0, v\dBR = 0, d(x, t) - 62 as |x| -> +00 (6.11) 

and modified initial conditions 

p(x,0) =po G 1^1'°°(Sh), 0<p<po<p<oo, (6.12) 

v(x,0)=VoGX„, d(x,0)-e2GH3(M2)^ (6.I3) 

where (v, ^tv) = in {R'^\Br) x / in flOOD . Here we have used the fact that v G L'^{I, H(^(5r)) 

is equivalent to v G L'^{I, H^(R2)) and v = in {R'^\Bn) x /. Hence one has v G C^{I, B.\R^) n 
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L°°(M2)) and dtv G L^{I, B.\R'^)) by zero extension, if v G C(J, X„) and ^t^f G L2(J,X„) with 

n = BR. 

Similar to the arguments in Section \AA\ we shall derive the uniform-in-time energy estimates 
in order to extend the local solution globally in time. This can be done by modifying the 
derivation of the uniform estimates in Section I4.4[ The only difference in arguments arises from 
the energy estimates for d in the whole space, recalling that we have shown the uniform estimates 
on II Ad||j;^2(7 l2(-q)) and || Ad||j;^oo(-7' l2(q)) in fl4.39p and fl4.46p in a bounded domain. In the case 
of r2 = M^, however, if d2 satisfies some geometric angle condition, then we can use the rigidity 
theorem to deduce the uniform estimates or i || A d|| r^ni.^mn] and || Ad||ioo(7 l2(ir2)). The rigidity 
theorem, which was recently established in 20|, reads as follows. 

Proposition 6.2. Let > 0, cq > 0. Then there exists a positive constant ujq = w{d2,co), such 
that the following holds. 

Ifd := (^1,^2) : ^2 §1, Vd G H1(R2) with ||Vd||L2(M2) < Cq and ^2 > ^2; then 

||Vd||^4m2) < (1 - ^o)\\^d\\l, 



Consequently, for such a map the associated harmonic energy is coercive, i.e.. 



|Ad+ iVdrdI 



> f^dlAdll^. 



+ llVdl 



(6.14) 



On the other hand, we have used the fact d = 1 and the boundedness of domain Q to get the 
uniform estimate ||d||j;^oo(/ l2(q)) in f l4.46p . For the case of the whole space, we need the uniform 
estimate ||d — e2||Loo(/ l2(]j2)), which can be deduced from d = 1 and (16. ip . More precisely, we 
have 



1 „ , 
— r d 



< 26 



e2|lL2 



+ 0||Vd| 



-I 
4 ' 



divv||^2^T[f2^ ] llVd"^ 



L2 



e2pdivv 



dx 



which, together with Gronwall's inequality, yields 



e2||L2 



< ||do - e2 



|2 

Il2 



02 

+ (4^T + -||divv||i.(,^^.(K.)) 



|Vd| 



L°°(/,L2 



(6.15) 



Consequently, plugging (I6.14p . (I6.15P and Proposition 16 . 1 1 with rfpg > into the deduction of 
the global estimates in Section 14. 4[ we immediately establish the following global existence of a 
unique solution (p„, v„,d„) to the approximate problem (I6.8p - (l6.13p : 

Proposition 6.3. Let Q = Br with R> 0, ^02 ^ positive constant, {P,e,p,p) satisfy ij.j'T^ 
and 5 G (0, 1]. Assume po satisfies (4^4^ > ^0 ^ X„ and 



do2 > ^02, |do(x)| = 1, do(x) - 62 G H2(r2). (6.16) 

Then there exists a unique triple {pn, v„, d„) defined on x / with the following properties: 
(1) {pni^n) satisfies the regularity lj^4.51\ ), and dn satisfies 

d„ - 62 G C°(/, H2(M2)) n L\I, H=^(R2)), dtd.r, G C°(/, L\R^)) n L2(/, h1(M2)). 
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(2) v„,d„) solves h6.9\) a.e. in Qt and l{6.10\) a.e. in x I, and satisfies lid. 8\) . 

(3) (p„,v„,d„) satisfies the finite and hounded energy inequalities as in Proposition \4-l\ with 

in place ofQ, where {pn,^n,'V Pn) = in (M^\r2) x /. Moreover, we have the uniform 
estimates ^-56^ , j-j-Sl^ , and 

|d„,| = 1 m X /, 

||V^d„||L2(]K2x7) + ||Vd„||i4(]g;2x/) < ^(^025 sup £^5(po, qo, do)). 

<56(0,1] 

||5tdn||L4/3(j^L2(m) + ||9td„||L2(j^(Hi(Q))*) <G{d^2^ sup £s{pQ, do) ■ 

<5e(o,i] 

As a consequence of Proposition 16.31 similarly to Section |5l we can use the standard three- 
level approximation scheme and the method of weak convergence based on Proposition 16.31 to 
establish the existence of weak solutions to the following problem: 

dtp + div(p v) = in X /, (6. 17) 

dt{p^) + div(pv (g) v) + VP(p) 

1, 



pAv + (p + A)Vdivv-z/div(^Vd0Vd-^|VdplJ in 5^? x J, (6.18) 
9td + v- Vd = e(Ad+|Vdpd) in x /. (6.19) 
More precisely, we have the following conclusion: 

Proposition 6.4. Let the initial data (po, mo) satisfy and ^1.10\) with Br in place of VL, 

and do satisfy 

do2 > do2, |do| = 1 a.e. m do(x) - ea G H^(R^). (6.20) 

Then the initial- boundary value problem Ili6.1'7\ )- ^6.19\) has a global weak solution (p, v, d) defined 
on X / for any given T > 0, such that 

(1) (p, v) satisfies the regularity U.13\) and (TJ^ with in place of Q, and d satisfies the 
regularity /ll.26\) . Moreover, (p, v) = m (M.'^\Br) x /. 

(2) Equations ^6.17\ ), ^6.18\) hold in V\Br x /), the equation ^6.19\) holds a.e. in R^ x /, the 
equation ( [g. i7| ) is satisfied in the sense of renormalized solutions, and the solution satisfies 
the energy equality U.20\) with in place offl. 

(3) Additional estimate: 

f f ^(|Ad|2 + |Vd|^)dxds < So, 
Jo Jr^ 2 

where the constant vjq depends on c?02) ^ ^'^^ 

^0=/ fi^^^l{p„>o} + (5(po)V^+ / ^^^5^dx. 
JBr \^ Po / Jr2 I 

Remark 6.1. To replace f l6.16p by (16. 20 p . we have used the fact that there exists a sequence of 
approximate functions {dJJ^}^^]^ C H-'^(R^) +62, such that 

|d^| = 1, d^ - do ^ in H1(R2), and c?^2 > ^02/2. 
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Finally, we can follow the arguments in 3J, Section 7.11] for the compressible Navier-Stokes 
equations to prove Theorem 11.21 by using Proposition 16.41 on the bounded invading domain 
and letting i? — )■ oo, where one should use the embedding theorem 



H^ip!) ^ L^ip!) for any p>\ and any bounded Lipschiz domain Vt' C M^, 

to replace L®(M'^) )■ D-'^'^(]R'^) in [34, Section 7.11] in the treatment of the velocity v. Here we 
omit the proof, since the additional limit process on d is trivial. We mention that, in the limit 
process on d, we also have the weak convergence as in (15.21) . and the strong convergence as in 
fl5.1l) and (15. 3p with any bounded space C M^. 
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